Dissipation and quantum phase transitions of a pair of Josephson junctions 
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A model system consisting of a mesoscopic superconducting grain coupled by Josephson junctions 
to two macroscopic superconducting electrodes is studied. We focus on the effects of ohmic dissi- 
pation caused by resistive shunts and superconducting-normal charge relaxation within the grain. 
As the temperature is lowered, the behavior crosses over from uncoupled Josephson junctions, sim- 
ilar to situations analyzed previously, to strongly interacting junctions. The crossover temperature 
is related to the energy-level spacing of the grain and is of the order of the inverse escape time 
from the grain. In the limit of zero temperature, the two-junction system exhibits five distinct 
quantum phases, including a novel superconducting state with localized Cooper pairs on the grain 
but phase coherence between the leads due to Cooper pair cotunneling processes. In contrast to a 
single junction, the transition from the fully superconducting to fully normal phases is found to be 
controlled by an intermediate coupling fixed point whose critical exponents vary continuously as the 
resistances are changed. The model is analyzed via two component sine-Gordon models and related 
Coulomb gases that provide effective low temperature descriptions in both the weak and the strong 
Josephson coupling limits. The complicated phase diagram is consistent with symmetries of the two 
component sine-Gordon models, which include weak to strong coupling duality and permutation 
triality. Experimental consequences of the results and potential implications for superconductor to 
normal transitions in thin wires and films are discussed briefly. 



I. INTRODUCTION 

Understanding the effects of dissipation on quantum 
phase transitions has proved to be a challenging problem 
in many contexts including quantum Hall transitions,^ 
and quantum critical points in antiferromagnets.^ Transi- 
tions from superconductor to "normal" metal or insulator 
in thin wires and films have been extensively studied, '^^^ 
as well as in Josephson junction arrays^"^^ and super- 
conducting nanowires.^'*"^* One of the most intriguing 
aspects of these transitions is the role of dissipation. ^^"^^ 
Theoretically, there has been extensive work on the ef- 
fects of dissipation on a single resistively shunted Joseph- 
son junction (RSJJ). The resistor can be modeled theo- 
retically as a Caldeira-Leggett ohmic heat bath,^^"^^ and 
precise predictions for the transport properties can be 
worked out (see [33] for a review). The system under- 
goes a superconductor-to-normal transition at zero tem- 
perature when the shunt resistance increases through a 
critical value equal to the quantum of resistance Rq — 
h/Ae^ — 6.53 fcfi. Recent experiments by Penttila et.al.^'^ 
showed good agreement with the theoretical analysis. 

Arrays of RSJJ have been studied in the same frame- 
work in terms of the local physics of the individual 
junctions. By percolation arguments, this local 
physics has been argued to apply to granular films and 
wires with the superconductor-to-normal transition in 
these extended systems occuring when the individual 
shunting resistances along a critical percolation path be- 
come equal to Rq.^^ 

The prediction for destruction of superconductivity 
via this local mechanism is in striking contrast to what 
one would expect in the absence of dissipation: domi- 
nation near to the quantum phase transition by collec- 
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FIG. 1: A mesoscopic superconducting grain connected to 
superconducting leads via Josephson junctions and resistive 
shunts. 



live long-wavelength quantum fiuctuations rather than 
local physics. In addition to the nature of the transi- 
tion, where it would be expected to occur as parame- 
ters of the system are varied is strikingly different for 
the two pictures. The long- wavelength quantum fluctu- 
ations should be controlled by the interplay between the 
Josephson couplings among grains and the Coulomb in- 
teractions, the former acting to decrease the phase fluctu- 
ations and the latter to decrease the charge fluctuations. 
If long-wavelength physics dominates, the location of the 
transition would thus be expected to depend markedly 
on the strength of the Josephson couplings. In contrast, 
for a single junction, and by naive extension for a net- 
work of junctions, the location of the dissipation induced 
transition would be entirely determined by the shunting 
resistances, independent of the Josephson couplings. 

The primary purpose of this paper is to begin to recon- 
cile these two approaches by studying a deceptively sim- 
ple system: two resistively shunted Josephson junctions 
coupled in series through a superconducting grain. This 
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system, in addition to its intrinsic interest, provides a 
simple paradigm for the competing effects of dissipation 
and quantum fluctuations on superconductivity. 

An important simplification in all previous theoretical 
studies of JJ arrays is the assumption that the supercon- 
ducting grains are sufficiently large that they can effec- 
tively be treated as macroscopic. In the case of several 
junctions in series, such an assumption leads to a result 
that the superconductor-to-normal transition occurs on 
each junction separately and takes place when the values 
of the individual shunting resistances are equal to the 
quantum of resistance Rq = h/{2e)'^. In this paper we 
take into account the effects of finite size grains, specifi- 
cally by considering two bulk superconducting leads con- 
nected by a pair of Josephson junctions in series through 
a mesoscopic grain. We show that the quantum dynamics 
of the two junctions, which are independent over a wide 
range of temperatures, become strongly coupled below a 
characteristic crossover temperature. In the low temper- 
ature regime, this simple system exhibits surprisingly rich 
behavior, including two distinct superconducting phases. 
In some regimes of parameter space, the superconductor- 
to-normal transition between the two macroscopic leads 
is determined by the total shunting resistance of the 
system, rather than individual resistances of the junc- 
tions; while in other regimes its location depends on the 
strengths of the Josephson couplings as well as the shunt- 
ing resistances. In this latter case, the corresponding 
critical behavior becomes very different from the single 
junction case. 

The basic system is shown in Fig. I. Dissipation oc- 
curs in ohmic shunts between the superconducting con- 
tacts and the grain. Such systems may be understood in 
terms of a two-fluid model in which Cooper pairs tunnel- 
ing across Josephson junctions represent the superfluid, 
and electrons flowing through the shunt resistors repre- 
sent the normal fluid. ^^"^"^ The presence of two fluids in 
the middle grain suggests considering it as a double grain 
with a superconducting part and a normal part as shown 
in Fig. 2. We assume for simplicity that the normal 
and superconducting charges of the two parts experience 
the same electrostatic potential as they overlap in space. 
The chemical potentials of the two parts, however, do 
not have to be the same. When these differ, the result- 
ing electrochemical potential difference can cause charge 
relaxation within the grain that will act to equilibrate 
its normal and superfluid components. In this paper we 
assume a simple ohmic model of this relaxation with the 
conversion current 



where Vn and Vs are the electrochemical potentials of 
the normal and the superconducting fluids on the grain. 
The coefficient r is a phenomenological parameter of 
our model that we will call the conversion resistance. 
Decoupling of the two chemical potentials is similar to 
the nonequlibrium state of the superconducting and nor- 
mal fluids, as discussed for phase slip centers at finite 
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FIG. 2: Effective circuit consisting of two Josephson-junctions 
{Jl, J-z) connecting the macroscopic electrodes {(t)i, <f>2) to a 
mesoscopic grain. The grain is modeled in a two-fluid man- 
ner, as a superconducting grain {(jig) connected through a 
phenomenological resistance r to a normal- fluid grain {■^). 
Ri, R2 are the shunt resistors connecting the normal-fluid 
of the grain to the superconducting contacts, in which the 
normal-superconducting relax;ation is fast. 

current. '''"'"^^ We assume that the two leads are macro- 
scopic, so that there is perfect coupling between the su- 
perconducting and normal fluids in each of them (this 
corresponds to the conversion resistances in the leads be- 
ing negligible). 

The model we arrive at using the arguments above is 
quite general. One could also obtain it by considering 
the electromagnetic modes that Cooper-pair tunneling 
events excite as discussed in Appendix B. This alterna- 
tive approach does not require a two-fluid picture. 

It is worth pointing out that our system bears some 
resemblance to the Cooper pair box systems studied re- 
cently in the context of quantum computing and meso- 
scopic qubits.*^"'''" The charge on the grain could be used 
as a the quantum number of a qubit. The biggest obsta- 
cle to quantum computation is then the limited lifetime 
of the quantum state of the qubit. Quantum fluctuations 
and interactions with the environment limit the life time 
of such a state, so practical realizations of qubits require 
systems with low dissipation. In this paper, in contrast, 
we study the Cooper pair box system in a highly dissi- 
pative environment. 

This paper is organized as follows. In Section II we 
present a microscopic Hamiltonian and derive the quan- 
tum action. To ascertain the consistency of this deriva- 
tion we demonstrate in Appendix A that the classical 
equations of motion obtained from the action correspond 
to the electrodynamics of the circuit in Fig. 2. From the 
analysis of the quantum model we show the existence of a 
new temperature scale T* set by the level spacing in the 
grain. At temperatures higher than T* the two junctions 
are decoupled and can be considered separately. If the 
grain is macroscopic, T* ^ and the system is always in 
the decoupled regime. This is the case considered in the 
literature thus far.^^~*° For temperatures below T* one 
cannot neglect interactions between the junctions, and 
the effective low temperature description is given by two 
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coupled quantum sine-Gordon models. 

In section III we use renormalization group (RG) meth- 
ods to analyze the two component sine-Gordon theory 
in the limit of weak Joscphson coupling and obtain its 
phase diagram. We show that the system can have five 
distinct phases: fully supercondicting, FSC, where both 
jimctions arc superconducting; normal, NOR, where both 
junctions are normal and there is no phase coherence be- 
tween the leads; N1-S2, where junction one is normal and 
junction two is superconducting; S1-N2, where junction 
one is superconducting and junction two is normal; and 
SC*, in which Cooper pairs are localized on the grain, so 
individual junctions are insulating, but there is supercon- 
ducting coherence between the leads due to cotunneling 
processes. We provide simple arguments for the phase 
boundaries based on electrical circuit considerations of 
the effective shunting resistances for various Cooper pair 
tunneling events. 

In Section IV we analyze the system in the oppo- 
site regime of strong Josephson couplings using a dual 
two component sine-Gordon model and considerations of 
quantum phase slips. The RG analysis is again supple- 
mented by effective shunting resistance arguments which 
determine the action of the various quantum phase slip 
processes. It is found that the phase diagrams obtained 
in the weak and strong coupling limits differ in the loca- 
tion of the NOR to FSC phase boundary. 

In Section V we show that the difference between 
strong and weak coupling phase diagrams signals the ex- 
istence of a novel regime with the fully normal to fully 
superconducting transition controlled by a critical fixed 
point at intermediate Josephson coupling. We analyze 
the appropriate fixed point, whose properties depend 
continuously on the resistances, and discuss the RG flow 
in its vicinity. 

In Section VI we explore the surprisingly rich symme- 
tries of the two-junction system. In addition to a weak- 
to-strong duality, the system also exhibits a permutation 
triality that implies that aspects of the phase diagram 
are invariant under interchange of any of the three resis- 
tances involved in the dissipative transport. 

In Section VII we review some experimental implica- 
tions of our work and discuss such questions as observa- 
tion of the crossover temperature scale T*, the experi- 
mental identification of the novel superconducting phase 
SC*, and the universality of the resistance at the super- 
conductor to normal transition. We also suggest that 
our results may be relevant for understanding some puz- 
zling experimental results on superconductor to normal 
transitions in thin wires and films. 

Finally, in Section VIII we summarize the main results. 
To maintain the coherence of the presentation we dele- 
gate most of the technical calculations to appendices. In 
particular, the renormalization group analysis of the two 
component sine-Gordon model and the relations to clas- 
sical Coulomb gasses are given in Appendices D (weak 
coupling), and E (strong coupling). 



II. MICROSCOPIC MODEL 
A. Hamiltonian of the two-junction system 

The system we wish to describe consists of a meso- 
scopic superconducting grain situated between two 
macroscopic superconducting leads (Fig. 2). The 
grain interacts with the leads both electrostatically and 
through a weak link. The electostaticl interaction is ca- 
pacitative while the weak link allows the flow of both 
Cooper-pairs and normal electrons. Cooper pairs flow 
through a Josephson junction from the superconducting 
part of the grain to the leads. Normal electrons flow from 
the normal part of the grain to the leads through what 
we model as a shunt resistor. 

In order to understand the quantum dynamics of this 
system we must first obtain an appropriate low-energy 
effective Hamiltonian. This should include the charging 
energy for the grain and leads, the Josephson coupling 
energies for the junctions, and appropriate Hamiltonians 
for the shunt resistors which can be approximated by 
heat baths. 

The charging energy of the system includes both elec- 
trostatic and electrochemical capacitances. All the is- 
lands (here we use the term island to denote either the 
electrodes or the grain) have part of their charge, Qsi, 
in the form of superconducting Cooper pairs and part 
of their charge, Qjvi, in the form of normal fluid. Both 
kinds of charge contribute to the electrostatic potential 
and have their own compressibility. The electrochemical 
potentials for the superconducting and normal electrons 
on island i are 

Vsi = + DsiQsi , ^ 

(2) 

Vivi = <pi -|- Dj^iQ^i- 

The index i is summed over electrodes 1,2, and the grain 
g; (pi is electric potential; -Di's are the inverse of the com- 
pressibilities of the fluids S and N, in a non-interacting 
approximation, e^DNi is the level spacings of the normal 
electrons in the island i.^^ The electrostatic potential on 
island i is related to the charges on all the islands via the 
capacitance matrix dj : 

Vi = 12C!r.\Qsj + QNj). (3) 

j 

Hence, for the electrochemical potentials we have 
Vsi = HM'slfisj + C^'QN,) 

(4) 

where we deflned 

>^mj = +DNiSij 

with Sij a Kronecker delta. By integrating out the 
electro-chemical potentials in (4) we find the charging 
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part of the Hamiltonian 

-1 

+ QsiQnj- 

At this point we introduce superconducting phases, (pi, 
on the islands and "normal phases" that we define 
formally to be conjugate to Qni'^^ 

[QNi,tp]] = -ieSij {Qsi,4>i] = -2ieSij 

(7) 

[QNi,(l>i]=0 [Qsi,iJj] = 0. 

By using (6) and (7), it is easy to verify that the Hcisen- 
berg equations of motion for the two phases give the cor- 
rect Josephson relations: 

The other important energies involving the supercon- 
ducting degrees of freedom are the Cooper pair tunnel- 
ings, with 

Hj = - Ji cos{(f)g - (pi) - J2 cos(^!>2 - (pg). (9) 

The dissipation in the ohmic shunts, Ri, R2, and the 
internal charge relaxation, r, are modeled following 
Caldeira and Leggett (see refs.^^'^^''^^ for a review). In 
this approach, the shunting resistances are replaced by 
collections of harmonic oscillators (heat baths), with ap- 
propriately chosen spectral functions: 

Hdis = nbath{Rl,2iJl - 2lPg) + Hbath{R2,2i}2 - 2iPg) 

+'Hbath{r, (pg -2ipg). 

(10) 

We will not give the explicit form of the appropriate 
Hamiltonians here, but in the next subsection we give 
the eff'ective actions obtained after integrating out the 
heat-bath degrees of freedom. The heat bath model is 
the simplest quantum model that gives the correct clas- 
sical equations of motion for systems with dissipation. 
Later in this paper we will discuss some of its drawbacks, 
however, we believe that it gives a qualitatively correct 
picture for a general mechanism of dissipation. 

Collecting all the terms, we obtain an effective Hamil- 
tonian that describes the system shown in Fig. 2: 

n{QNU QS^, <P^. ^r) =nQ+Hj+ H^s- (n) 



B. Imaginary time action 

From the Hamiltonian (11) and commutation relations 
(7), we can construct the imaginary time action and par- 
tition function for the system in Fig. 2: 

Z = JVQmVQsiV(PiVi)iexp{-S) 

S = -TeE^Io dr Qsi 0i - i Ei Jo dr Qm ^^ (12) 

-t- Jq dTn{QN2,Qsi,4>ii'4'i)- 

It is important to point out that in the presence of ohmic 
dissipation the phase variables (pi and tpi should be peri- 
odic at T = and t = i3 with no phase twists by multiples 
of 2tt allowed. This follows from the fact that a 2tt phase 
twist causes dissipation and is thus measurable. The 
ohmic dissipation allows continuous charge transfer (as 
opposed to transfer of multiples of e) from the shunting 
resistors to the grain. Therefore any non-integer charge 
induced by the gate voltage can be screened out. (For a 
more detailed discussion see^^'^^). This potential draw- 
back of the Caldeira-Leggett model of dissipation may 
be overcome if one introduces a more complicated form 
of dissipation, such as via quasiparticle tunneling (see 
e.g.^^). 

The quantum action in (12) is quadratic in Qsi and 
Qni, so they may be integrated out (for details, see 
Appendix Al). The electrochemical contribution is (in 
terms of the electrochemical potentials) 

= lo dT whw (Si (^Qii^Si - Vmf 

. (13) 

+ T,ij{siVsi + ViVNi) (sjVsj + VjVNj)j ■ 

This is very easy to interpret. The level spacings give rise 
to the first term in the brackets making a potential dif- 
ference between the two fluids on one island energetically 
costly. The second term in the brackets is the charging 
energy one would expect from a conventional system of 
islands, but the potential on each island is replaced by 
a weighted average of the normal-fluid potential and the 
superfluid potential: Vi = SiVsi + rjiVNi- 

In terms of the phase variables, the full action can be 
written as 
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Z = J V(j}iVil)i exp {-Sq -Sj- Sdis) 

^7 = dT{-Ji C0s{4>g - (j)i) - J2 C0s((/)2 - (l>g)) 

Sdis = /3E.„ ^ (W|2V'i,(a,„)-2V',,(.„)|^ + ^|2V'2,(.„)-2V',,(a,j|^ + ^|</',,(.„)-2V',,(.„)|^) 



where the Matsubara frequencies are uin = 27rTn, and we 
have defined 

CQi = {Dsi + DNi)-' (15) 
Si = DNi/{Dsi + DNi) (16) 

and 

% = Dsi/{Dsi + Dm) ■ (17) 

An important consequence of the domain of the phase 
fields (j}g and ijj being the real line rather than a circle, is 
that the Berry phase has no ofFcct on the behavior of the 
system. A Berry phase could arise if we included the gate 
voltage effects in (2) and (6) by shifting Qsg Qsg — Qo, 
which would lead to additional terms in the action (13) 
of the form iQq <j)g . But because 27r phase twists are 
not allowed, the additional action vanishes due to the 
periodic boundary conditions in imaginary time. 

As a consistency check on the action (14), we demon- 
strate in Appendix A 2 that its real time equivalent gives 
rise to equations of motion that coincide exactly with the 
basic electrodynamic equations for the circuit in Figure 
2. 

In this paper we consider the limit of macroscopic elec- 
trodes, so we can set the corresponding _Di = D2 = on 
these. The first term in (13) then imposes perfect cou- 
pling between the superconducting and normal fluids in 
the electrodes, i.e., (pi = 2^)1 and <j)2 = 2^2- Note that 
this assumption does not restrict us to taking an infi- 
nite capacitance for the electrodes: the inverse of the 
level spacing grows as the volume of the grains, whereas 
capacitances increase only linearly with the dimensions. 



We restrict our discussion to the case when the largest 
capacitances in the system are the mutual capacitances 
between the electrodes and the grain, Ci and C2, for 
electrodes one and two respectively. In Appendix A 3 we 
show that in this case the charging part can be simplified 
if we introduce the phase difference variables 

Al = 0g - (pi 

A2 - 02 - 09 (18) 

Ag = (j)g- 2^g, 

and the center of mass variable, 

^ ^ Cll-f-Cl2+Cl, , ^ C22 + Cl2 + C2g , 

I Clg + C2g+Cgg „ J, I Clg + C2g+Cgg n I 

Ctot 9 ^9 Ctct '19 ^''^9' 

where 

Ctot = Y,Cij (20) 

(note that Ctot is not affected by the mutual capacitances 
Ci and C2 but is determined by the capacitance of the 
system to the ground). We thus have 

+C2iA2 + VgKf + CQAl + Ctot^^) (21) 

The center of mass coordinate, completely decouples 
from the phase differences in the charging part of the 
action, and it is not present in Sj and S^ns ; these can be 
written as 



Sj = dT{- JiCOs{Ai) - J2C0S{A2)) 



(22) 



Therefore the center of mass coordinate, factors out in the partition function. From Eqs. (21)-(22) we see that Ag 
appears quadratically in the action and can be integrated out. After this integration, and also after neglecting terms 
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involving Ci/Cq, C^/Cq <C 1, we obtain 



Looking at the ubiquitous denominators of Eq. (23) we 
notice the expression: 



h n 1 1 



The scale for the Matsubara frequencies, w„, is set by 
temperature, hence a new temperature scale emerges 
from (23): 

= {2efiDs + Dr,)RQ (^1 + i- + i-^ . (24) 

This is the level spacing on the grain (1/Cq = Ds + D^) 
times a dimensionless resistance dependent factor, and it 
is also of the order of the inverse escape time from the 
grain. 

High temperature limit. When T ^T* the denomina- 
tor in (23) is dominated by \uj\ ^ T*, and the effective 
action for high temperatures is 

5 - ^/?E (i I Ai, (.„) I' + C,u;l/h) + 



h 1^2, (a-„) 



Al,(c.„)A2,(- 



(25) 



lo'l / h/CQ(l + \u:„\rigCiRi/h)(l + \u„\rigC2R2/h) \\ 

R1R2 \ R J J 

+Sj. 



In this limit we see that the interaction term between the 

two junctions (which is T-independent to leading order 
in Ci/Cq, C2/CQ) is negligible compared to the other 
resistive and capacitative parts of the action; the two 
junctions are thus effectively decoupled for T ^ T*. The 
dissipations for the two junctions in this limit are set 
simply by the individual shunt resistances Ri and i?2- 
This is the limit that has been discussed in the litera- 
ture; its validity at low temperatures relies on the basic 
assumption of m,acroscopic grains, for which T* = 0. 

Low temperature limit. At temperatures T below T* 
(we assume that T* < n/{RiCi) and T* < h/{R2C2)) a 
qualitatively different picture emerges in which coupling 
between the two junctions becomes important. The low- 



energy effective theory is 



Sj = 



Sc = PE 

Sd 



dT{-Ji cos(Ai) - J2 cos(A2) 
-J+cos(Ai + A2)) 

2 



2(2e)2 

^ At G A 



ri?.2Ai+rJ?.iA2 
rRi+rR2+RiR2 



(26) 



with 



and the matrix 



G = 



A= (Ai,A2) 



Rq_ fr + R2 r 
2nY \ r r + Ri 



where 



Y = rRi + rR2 + R1R2. 



(27) 



(28) 



(29) 



In the equations above wc have added, for future pur- 
poses, a lead-to-lead Josephson coupling term represent- 
ing cotunneling processes via the grain. This term de- 
scribes a Cooper pair tunneling (pair-tunnel event) from 
the left electrode to the right electrode (see Fig. 8) via a 
virtual intermediate state with an additional pair on the 
grain. Such processes appear perturbatively at second 
order in Ji and J2 and will be generated in the RG flows 
for the action (26) (see discussion below Eq. (30)). 

It is important to note that level spacing Dg = Dgg + 
Dng only appears in Sc via the quantum capacitance 
Cq = Dg^: whose precise form will not matter except 
to yield a high frequency cut-off. By the same token, a 
different form of the capacitative energy of the leads and 
the grains would only modify Sc and not change any of 
the analysis presented in this paper. 

Action (26) is one of the main results of this paper, and 
in the following sections we will mostly be concerned with 
studying its properties. 
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Superconducting Grain 



FIG. 3: Physical interpretation of expanding the Ji cos Ai 
terms in the action (26). The weak Joscphson coupling action 
can be mapped to a theory of interacting Cooper-pair tun- 
neling events (Coulomb-gas representation) with each pair- 
tunnel "charge" corresponding to a Cooper-pair transferred 
through one of the junctions. The co-tunneling events which 
transfer Cooper-pairs from lead to lead are also shown; in the 
Coulomb-gas representation these correspond to pair-tunnel 
dipoles. 

III. WEAK COUPLING ANALYSIS 

In this section we analyze the low energy properties of 
the system in the weak Josephson coupling limit. 

A. Renormalization Group equations 

In the limit of weak Josephson couplings {Ji}, the 
quantum action (26) can be analyzed directly in the gen- 
eralized sine-Gordon representation. In Appendix D2, 
RG flow equations are derived to second order in the 
Josephson couplings: 

^ = J2(l-^) + ^JiJ+ (30) 

In writing these we have set a combination of the short- 
time cut-ofFs to be equal to one. In physical imits, the 
energy cut-off is of order the charge relaxation rate of 
the junctions in units of which we are here measuring 
the {Ji}. The first order terms in the RG flows arise, as 
usual, from integrating out fast modes in the quadratic 
part of action (26). The second order terms are obtained 
from recombinant terms in the expansion in powers of J's 
of (26). These can be understood physically; pair-tunnel 
events on junctions one and two can combine to form a 
cotunncling event between the two leads, while a cotun- 
neling event plus a pair-tunnel in the opposite direction 
across one of the junctions is equivalent to pair-tunneling 
across the other junction (for details see Appendix D2). 
Prom (30) we see that, as claimed in the previous section, 
J+ gets generated at low energies even if we start with a 
model in which = 0. 



B. Weak coupling phase diagram 

Surprisingly, the simple flow equations (30) give rise to 
five different regimes. When all J's are irrelevant about 
the uncoupled fixed line so that they flow to zero, the 
system is in the normal state with no supercurrents be- 
tween the leads or between either lead and the grain. 
This normal (NOR) phase occurs if Ri and R2 are both 
sufficiently large. When all J's are relevant and grow 
under the RG flows, the systems is in a fully supercon- 
ducting phase that we denote FSC. This occurs if all the 
resistances are sufficiently small. For intermediate ranges 
of the resistances, the situation is more complicated. 

When only one out of the three Js is relevant while 
the other two flow to zero at low energies, the system is 
in a "mixed phase" ; as we shall see, there are three such 
phases. When the only relevant coupling is Ji, junction 
one is superconducting and junction two is normal, we 
call this phase S1-N2. With respect to lead-to-lead trans- 
port this is like the normal phase. Analogously we will 
have an N1-S2 phase when J2 is relevant but Ji and J+ 
are not. Rather surprisingly, there can also be a situation 
in which J+ is relevant but Ji and J2 are not. This is a 
phase in which individual junctions are normal, but the 
circuit as a whole is superconducting and Cooper pairs 
can flow freely between the leads. We denote this phase 
SC* . Physically, it corresponds to Cooper pairs being lo- 
calized on the grain, so that the individual junctions are 
normal, however the cotunneling processes, via virtual 
Cooper pair excitations on the grain, induce supercon- 
ducting coherence between the leads. A similar phase 
was discussed by Korshunov^^ and Bobbert et. al.^° in 
the context of one dimensional Josephson junction arrays. 

Inspection of the flow equations shows that as long as 
Ri,R2,r > there cannot be phases in which two of the 
three Js grow while the third flows to zero: the coupling 
terms in (30) from the two growing ones will drive the 
third J to grow as well. The system will then be in the 
fully superconducting (FSC) phase. 

To lowest order for small Js, the phase boundaries be- 
tween these phases are set by the relevance of Ji, J2, 
and J+ about the decoupled (normal) fixed line; these 
are determined by the combinations Ri +r, R2 + r, and 
Ri + R2 respectively, leading to the weak-coupling ap- 
proximate phase diagram shown in Fig. 4. 

This simple analysis, however, is not sufficient to ob- 
tain the correct phase diagram. For example, by looking 
at Fig. 4(a) we see that it has a superconductor to nor- 
mal transition of junction one inside the superconducting 
phase of jimction 2 (the N1-S2 to FSC transition). In this 
situation the RG equations derived in the vicinity of the 
uncoupled (normal) fixed line are no longer valid. 

A better approximation for the N1-S2 to FSC transi- 
tion can be obtained by noting that the fiuctuations in 
phase difference across junction two, A2, will be small in 
the N1-S2 phase. Thus in this regime we can approxi- 
mately set A2 = in (26). This modifies the RG fiow for 
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FIG. 4: Naive weak Josephson coupling phase diagram from first order RG. (a) For r < 0.5Rq four phases are present. In the 

superconducting (FSC) phase both junctions are superconducting, in the normal phase (NOR) both junctions arc normal. In 
the two remaining phases one junction is conducting and the other normal (N1-S2 and S1-N2). (b) When 0.5 < r/Rq < 1 there 
is a new phase, SC* , in which charges are localized on the grain. Both junctions arc thus normal, nevertheless Cooper-pairs can 
coherently tunnel between the leads (co-tunneling) so that this phase is superconducting with respect to lead-to-lead transport, 
(c) For r > Rq the FSC phase and the N1-S2 and S1-N2 phases disappear. The system is either in the SC* phase, or in the 
NOR phase. 



Ji to 




Ri 



rB.2 
r+R2 



Rq 



(31) 



[Wc will sec later that in the Coulomb gas language, equa- 
tion (31) corresponds to including the screening effects of 
unbound type two charges when considering the unbind- 
ing transition for charges of type one (see Appendix D 2 
for details).] From (31) we find that the N1-S2 to FSC 
boundary gets shifted to 



Ri + 



ri?2 



= Rq- 



(32) 



r + i?2 

Similar modifications of the phase boundaries appear for 



all transitions that involve ordering of one field in the 
presence of order in another: S1-N2 to FSC (ordering of 
A2 when Ai is ordered); and SC* to FSC (ordering of 
Ai and A2 when Ai -|- A2 is ordered). 



The corrected — and rather complicated — weak- 
coupling phase diagram is shown in Fig. 5. A particu- 
larly interesting regime occurs for r > Rq . In this regime 
the two-junctions cease to behave as such, instead, they 
behave much like a single junction shunted by the total 
resistance, Ri + R2, which therefore determines the lo- 
cation of the superconducting-to-normal transition be- 
tween the two leads. This will be discussed further in 
Section VII. 
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FIG. 5: Complete weak Josephson coupling phase diagram. Phase boundary formulas apply everywhere, although they are 

each given in only one graph, (a) When r < 0.57?q four of the five phases arc present; each junction is either normal or 
superconducting, (b) For larger r, the shape of the phase boundary between the FSC and NOR phases changes, (c) When 
0.75 < t/Rq < 1 the SC* appears and all five phases are present, (d) When r > Rq only SC* survives of the mixed phases, 
and FSC disappears. 



C. Circuit theory for weak coupling 

In this subsection we show how the phase diagram of 
Fig. 5 can be obtained by simple physical arguments. 
Before proceeding it is useful to recall such an argument 
for a single junction. 

Wc want to investigate the stability of the supercon- 
ducting state of a single Josephson junction with an 
ohmic shunt. In the superconducting phase, Cooper pairs 
are delocalized between the leads. Each Cooper pair tun- 
neling event changes the charge on the junction by 2e. 
This charge needs to be screened by the normal elec- 
trons in the shunt thereby causing a voltage drop to ap- 
pear across the junction. By the Josephson-relation, this 
voltage drop induces a change in the phase difference 
across the junction. The superconducting phase with de- 
localized Cooper pairs will survive only when the phase 
change due to one Cooper pair tunneling event is less 
than 27r (otherwise the phase becomes delocalized) . From 



circuit equations and the Josephson relation we find 

2e = / Iwdt = / -z— = — =- / -r-rfi 
J ^ J Rs 2eRsJ dt 



2eRs 



-Ac6 
'(33) 

where In is the normal screening current, AV is the volt- 
age difference across the junction, and Acj) is the phase 
change due to a Cooper tunneling. Rewriting the last 
relation as 



27r Rq 



(34) 



we obtain the usual condition; the shunted Josephson 
junction is superconducting when Rs < Rq. We can 
summarize this argument by saying that a Cooper pair 
tunneling event provides a current source with a mag- 
nitude that depends on the shunting resistance. By 
the Josephson relation, this leads to a phase fluctua- 
tion across the junction, and the superconducting phase 
is only stable when this phase fluctuation is less than 
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27r.[Notc that this argument docs not rcaUy yield the ex- 
act condition: a muhiphcative factor of order unity could 
have arisen. A fuller analysis, as from the RG flows, is 
needed to obtain the correct coefficient.] 

Applying this approach to the two junction system of 
Fig. 2 effectively reduces the problem to determining 
the effective shunting resistance associated with Cooper 
pair tunneling events in various situations. As in the sin- 
gle junction case, a Cooper-pair tunneling can be simply 
modeled as a current source. 

• To find the transition between S1-N2 and the NOR 
phase, consider a Cooper pair tunneling across junc- 
tion one with junction two insulating and acting as a 
circuit disconnect. The effective resistance that makes 
a circuit with the current source is then Ri + r (see 
Fig. 6(a)), and the phase boimdary is at i?i 4- r = Rq. 
Analogously for the N1-S2 to NOR transition we have 
the circuit shown in Fig. 6(b) and a phase boundary 
at i?2 + r = Rq. 

• The SC* to NOR transition is marked by the prolifer- 
ation of cotunneling processes in which a Cooper pair 
moves between the leads, but with both junctions in- 
dividually insulating. The circuit describing this case 
is depicted in Fig. 6(c), with the cotunneling process 
described as two current sources forcing the same cur- 
rent through both Josephson junctions. The cotun- 
neling process leaves no charge on the grain, and is 
hence screened only by normal currents flowing in the 
resistors Ri and i?2- The effective shunting resistance 
in this case is Ri + R2 and the phase boundary is at 
Ri + R2 = Rq. 

• The transition between N1-S2 and FSC occurs while 
junction two is already superconducting and can hence 
be replaced by a short in the circuit (see Fig. 7(a)). 
The effective shunting resistance across junction one 
then involves r and R2 in parallel, as well as Ri ; there- 
fore the phase boundary for this transition occurs at 
Ri + ri?2/(r + R2) = Rq- By the same token the 
transition between S1-N2 and FSC takes place when 
R2 + rRi/{r + Ri) = RQ. 

• To understand the FSC to SC* transition we need to 
consider the regime in which cotunneling maintains 
coherence between the leads; therefore these are ef- 
fectively connected by a short in the circuit as shown 
in Fig. 7(b). Now consider a Cooper pair tunneling 
from the grain to one of the leads, say two. The ef- 
fective resistance seen by a tunneling Cooper pair is 
r -|- RiR2/{Ri -\- R2) and the phase boundary is hence 
at r RiR2/{Ri + R2) = Rq- The effective shunt 
resistance for tunneling from the grain to lead one is 
the same. The nature of the SC* phase is as follows: 
Cooper pair tunneling events scramble phases across 
junctions one and two too much for the junctions to 
be coherent, so Cooper pairs become localized on the 
grain. Nevertheless, cotunneling events allow Cooper 
pairs to move between the leads, so there is a well 



, — WWV- 



I — 7\AAA^-^^AAMA — , 



I — ^^AAA 7WW — I 



FIG. 6: Effective circuits for pair tunneling events. A pair- 
tunnel corresponds to a current source, whereas a junction 
without pair tunneling acts as an open-circuit, (a) Effective 
circuit for a pair tunneling through junction 1. (b) Effective 

circuit for a pair tunneling through junction 2. (c) Effective 
circuit for a coherent load-to-lcad pair tunneling event. Since 
the current through junction 1 and 2 is the same, the resis- 
tance r is effectively disconnected in this case. 



defined phase difference between them that acts as a 
short between the two macroscopic leads as far as dis- 
sipation across the individual junctions. 

• The FSC to NOR transition line is, naively, a contin- 
uation of the S1-N2 to NOR and N1-S2 to NOR lines. 
This suggests that when considering fluctuations of the 
phase difference across junction one, we assume junc- 
tion two to be insulating, and vice versa. The consis- 
tency of such an approximation is highly questionable 
and reflects the limit of small J's as our starting point: 
by a weak-coupling analysis: for a weak-coupling limit 
to be valid, we should only approach phase boundaries 
from normal phases of the junction under considera- 
tion. 

It is worth pointing out that in all cases described 
above, the effective dissipation is decreased relative to 
that in the high temperature action (25). The most ex- 
treme case happens for the SC* to NOR transition which 
is determined by the total shunting resistance at low tem- 
peratures rather than individual resistances Ri and R2, 
which would determine the transitions between macro- 
scopic grains. In the SC* phase the whole system be- 
haves as a single junction, and the dissipation is deter- 
mined by the resistance across the whole of the system. 



IV. STRONG COUPLING ANALYSIS 

We have seen that much can be concluded from the 
weak Josephson coupling analysis, in particular the na- 
ture of the five possible phases and some of the transitions 
between them. Yet some of the transitions could only be 
understood via a hybrid analysis involving some large 
and some small couplings, and, as pointed out above, 
the FSC to NOR transition cannot be analysed in a con- 
trolled manner from a weak- coupling analysis. Even to 
solidify the identification of all of the superconducting 
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a. 



b. 





FIG. 7: (a) Effective circuit for a pair tunneling through 
junction-1 when junction-2 is superconducting (J2 is relevant 
about the weak coupling limit), (b) Effective circuit for a pair 

tunneling event through junction-1 (or 2) when J+ is relevant 
and coherent lead-to-lead pair tunneling events proliferate. 

phases, we really need to go beyond weak coupling: as 
soon as one or more of the Js grows without bound, the 
system flows out of the regime of validity of the RG flow 
equations used thus far and we must ask where it flows 
to. 

In this Section we turn to the limit of large Joseph- 
son coupling and attempt to analyse the phases, phase 
diagram, and transitions in that limit. 



of the superconductor to normal transitions can be de- 
scribed as binding-unbinding transitions of this two com- 
ponent plasma. 

In the strong coupling case we can write a Coulomb 
gas representation for the quantum phase slips instead of 
the Cooper pair tunneling events. The phase slips also 
behave as a two-component gas — phase-slips on the 
two junctions — with logarithmic interactions between 
them. When the phase-slips across a junction prolifer- 
ate, it becomes normal; if instead their fugacity tends to 
zero at low-energy-scales, the junction is superconduct- 
ing. Mathematically, the strong coupling case can be 
analyzed by performing a Villain transformation to rep- 
resent the partition function (26) in terms of two types of 
interacting phase slips. This classical Coulomb gas can 
then be transformed into a new sine- Gordon model that 
is dual to (26). Appendix C describes the details of such 
transformations. We find 

Z = I D[ei]jD[92]exp{-S) 
with 

s = /3j:\un\eZ„Me^^ 

- /(f dr (Ci cos(^i) + C2 008(62) + C- cos (^1 - ^2)) , 

(35) 

where 6 = {61,62) and 

- 1 1 ( r + Ri -r \ 
M = G-^ = — — (36) 

is the scaled resistance matrix. The variables (i, (2, C- 
are the fugacities corresponding to the three types of 
phase-slips: Ci across junction 1; (^2 across junction 2; 
and C- a combination of these that corresponds to a 
phase slip across 1 and a simultaneous anti-phase slip 
across 2, thereby slipping the phase on the grain with 
respect to both of the superconducting leads. 



A. Sine-Gordon action for quantum phcise slips 

When the Josephson couplings are large, the system is 
usually in the vicinity of one of the classical minima of the 
Josephson potentials so that Ai « 27rni, A2 « 27rn2 with 
ni and n2 integers. Only rarely does the system undergo 
a tunneling event in which one or both of the phases 
winds by 2tt. Such phase tunneling processes between 
minima of the classical potential are quantum phase slips 
(QPS).^^ When QPS across it are suppressed at low tem- 
peratures a Josephson junction is superconducting, but 
when they proliferate the junction is incapable of sup- 
porting supercurrcnts and becomes normal. 

In weak coupling wc analyzed the low energy action in 
terms of Cooper pair timneling events. As discussed in 
Appendix B, this is equivalent to a classical Coulomb gas 
with two types of charges corresponding to pair tunneling 
events through the two junctions; dissipation gives rise to 
effective logarithmic interactions among these. Various 



B. Phase diagram 

Following the steps leading to Eq. (30) we readily ob- 
tain the flow equations for the phase slip fugacities Ci , C2 
and 

*=^>-5%)^*«- (37, 

1f = <-(i-;Tt5) + K.&. 

where we use Y = R1R2 + rR\ + ri?2- These flow equa- 
tions are correct to second order in the Qs, being simply 
the analog of Eqs. (30) for the weak coupling limit. We 
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again work in units in which the short time cut-ofF — here 
related to the "transit time" for a least-action phase slip 
— is unity. 

Growth under renormalization of a fugacity Q corre- 
sponds to proliferation of the corresponding QPS and 
hence destruction of superconductivity across the respec- 
tive junction in the case of Ci or 1^2, or between the grain 
and the rest of the system in the case of C- • 

Equation (37) gives rise, as did the weak coupling anal- 
ysis, to five phases. When all C's are irrelevant and flow 
to zero, the system is in the fully superconducting state 
(FSC) since isolated phase slips all cost infinite action. 
Conversely, if all C's are relevant, we expect the normal 
state (NOR) to obtain. As in the weak coupling case, 
three mixed phases appear when only one of the fugac- 
ities is relevant. When Ci is relevant and ^2 and are 
not, the system is in the N1-S2 phase; analogously a rel- 
evant ^2 and irrelevant Ci and C- signal the S1-N2 phase. 

If ^_ is relevant but Ci and C2 are not, the special 
SC* phase occurs. In this phase only QPS dipoles pro- 
liferate; these consist of a phase slip across one junction 
and an anti-pha.se slip across the other. Isolated phase 
slips across individual junctions will not occur in the 
SC* phase. Superconducting phase coherence between 
the two leads is thus maintained, since the phase differ- 
ence between them is the sum of the phase differences for 
the two junctions, and a phase-slip on junction one gets 
canceled by its accompanying anti-phase-slip on junction 
two. But the phase difference between the leads and the 
grain is ill defined in SC* as a result of the proliferated 
QPS dipoles. We thus see that proliferation of the QPS 
dipoles induces charge localization on the grain. 

The transition between the two superconducting 
phases, SC* and FSC, is, from the point of view of phase 
slips on the individual junctions, a transition between a 
dipole-free state, FSC, in which all the phase slips will 
be bound in quadrapoles, and a phase, SC*, in which 
dipoles proliferate but single quantum phase slips still do 
not occur. Because of the free dipoles in the SC* phase, 
a single quantum phase slip between the two leads can 
consist of any combination of phase slips across the two 
junctions that add up to a total phase difference between 
the leads of 2tt. 

As in the weak coupling limit (Sec. IIIB), we could 
attempt to construct a naive phase diagram showing all 
five phases from the first order strong coupling flow equa- 
tions. This approach would give phase boundaries that 
depend onRi + ^,R2 + -^andr + (see Fig. 

8). But such an analysis, as in the weak coupling limit, 
is not sufficient: when one type of phase slip proliferates 
it will partially screen the interactions between the other 
types of phase slips. 

To do better we must consider the effects of the rele- 
vance of a C cos 9 term: this will cause the dual phase, 6, 
associated with the proliferating phase slips to become 
localized at an integer multiple of 27r. As the 6 will then 
not fluctuate appreciably about this at low energies, we 
can set it to zero. As for weak coupling, this suppression 



of some of the fluctuations will change the flows of the 
remaining fugacities and thereby modify the phase dia- 
gram. The complete phase diagram from such a strong 
coupling analysis is shown in Fig. 9. 

C. Circuit theory for strong coupling 

The strong coupling phase diagram of Fig. 9 can be 
simply interpreted in terms of the effective electronic cir- 
cuits. Although these arguments are dual to the ones 
used for weak coupling, we present them here for com- 
pleteness. 

Again it is useful to start by considering the case of a 
single junction, now starting from the superconducting 
regime. The normal state occurs when quantum phase 
slips proliferate. When a QPS occurs, the phase differ- 
ence across the junction changes by 2tt. By the Josephson 
relation, this generates both a voltage drop and charge 
flow through the normal shunt. In the normal state the 
Cooper pairs should be localized, therefore such a state 
can only be stable if the charge fluctuation caused by 
an individual QPS is less than 2e (again, the justifica- 
tion of the factor being exactly two really needing the 
fuller analysis). Prom Kirkhoff's laws and the Josephson 
relation we have 

2n= [ ^dt=^ [ Vdt = ^R, I Idt = ^Aq (38) 
J at h J n J h 

Here Ag is the amount of charge that passes through 
the shunt resistor as a result of the QPS. In units of the 
charge of a Cooper pair, 2e, this is 



We thus guess that the normal state is stable when Rg > 
Rq. The basic physics is that fluctuating QPS's act as 
voltage noise that gives rise to charge fluctuations on the 
junction. The insulating state is only stable when these 
charge fluctuations are sufficiently small: less than 2e. 

The generalization of the single-junction argument to 
the system in Fig. 2 requires analysis of the effective 
shunting resistances for the various QPS conflgurations. 
The quantum phase slips are effectively voltage sources. 
The phase slip dipole corresponding to is thus equiv- 
alent to two equal but opposite voltage sources across 
the two junctions so that there is no voltage between the 
two leads, but the grain is at a different voltage than the 
leads. 

• The FSC to N1-S2 transition is determined by the cir- 
cuit in Fig. 10(a). In this case junction two can be 
replaced by a short as it is superconducting on both 

sides of the transition. This gives an effective shunt- 
ing resistance i?i -|- ri?2/(J?2 + r) for the phase slip, 
so the transition occurs at i?i -|- ri?2/(J?2 + r) = Rq. 
Similarly, the transition between S1-N2 and FSC is 
determined by the circuit in Fig. 10(b), with the 
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FIG. 8: Naive strong Josephson coupling phase diagram from first-order RG flows. Phase boundary formulas apply everywhere 
although they are each given in one graph only, (a) In the range < r/Rq < 2/3, every junction is either normal or 
superconducting; (b) and (c) show the range 2/3 < r/Rq < 1 for which the SC* phase takes up some of the FSC parameter 
space, and also pushes the NOR phase towards the axis at the expense of N1-S2 and S1-N2. The dashed line is where the 
boundary of the insulating phase would have been if phase-slip dipoles had not been taken into account; (d) obtains in the 
range r > Rq for which the SC, N1-S2 and S1-N2 phases no longer occur, so that only SC* and NOR phases survive. 



effective shunting resistance at the transition being 
R2 + rRi/{Ri + r) = RQ. 

• To understand the FSC to SC* transition we need 

to consider a dipolc consisting of a QPS on junction 
one and a simultaneous anti-QPS on junction two, 
corresponding to a 27r phase twist on the interven- 
ing grain. In particular, wc need to know how much 
charge flows from the super electrons on the grain to 
the normal electrons on the grain during such a phase 
twist. An equivalent circuit is shown in Fig. 10(c), and 
we conclude that the phase boundary should occur at 
r -|- i?ii?2/(-Ri + ^2) = Rq^ as the charge must flow 
through r and either i?i or R2. 

• The transition between S1-N2 and the NOR phase is 
determined by the relevance of the QPS on junction 
one when junction two is insulating. The correspond- 
ing circuit is shown in Fig. 11(a); since the effective 
shunting resistance is Ri + r, we find a phase bound- 
ary ai Ri + r = Rq. Similarly, the N1-S2 to NOR 
transition is at i?2 + J' = Rq- 

• The transition between SC* and NOR is determined 
by the effective circuit in Fig. 11(b). In the SC* phase 



the component that is incoherent with the rest of the 
system is the grain. Since phase coherence between 
the leads is maintained, charge can flow freely from 
lead to lead via virtual super-conducting electrons on 
the grain unhindered by the phase fluctuations on the 
grain. But if some charge flows through r to the nor- 
mal electrons on the grain, this current will couple to 
the phase-slip dipoles and induce a large voltage drop; 
hence r becomes effectively a disconnect in the SC* 
phase. The destruction of lead-to-lead superconduc- 
tivity that characterizes the SC* to NOR transition 
thus occurs at Rx+ R2 = Rq. 



• The FSC to NOR line is naively a continuation of the 
S1-N2 and N1-S2 lines. This suggests that to approach 

the transition line from the superconducting side, when 
we consider a QPS in junction one we assume junction 
two to be superconducting and vice versa. This highly 
questionable approximation reflects the limitations of 
our strong coupling analysis for the FSC to NOR tran- 
sition; we will analyze it more carefully below. 
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FIG. 9: Complete strong Josephson coupling phase diagram. Phase boundary formulas apply everywhere, although they are 
each given in only one graph, (a) r = for which the two junctions are effectively independent; (b) range < v/Rq < 2/3; 
(c) range 2/3 < r/Rg < 3/4; (d) range 3/4 < r/Rq < 1 where all five phases are present (e) range r > Rq for which the two 
junctions act like a single junction with a shunt resistor Ri + i?2. 



V. INTERMEDIATE COUPLING FIXED POINT 



In the previous two sections we have analysed the zero 
temperature states and transitions between them in both 



the weak and the strong Josephson coupling limits. In 
both cases, we found that there were some regimes that 
could not be adequately analysed. In this section an- 
alyze the intermediate coupling behavior, finding that 
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FIG. 10: Effective circuits for transitions to the FSC phase 
in the phase slip picture. Phase sUps correspond to a voltage 
source across the corresponding junction, (a) phase slip on 
junction 1; (b) phase slip on junction 2; (c) slip-anti-slip pair 
which corresponds to slipping the phase of the grain relative 
to both the leads. 
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FIG. 11: Effective circuits for transitions to the NOR phase in 
the phase slip picture, (a) phase slip on junction 1 when junc- 
tion 2 is insulating (^2 is relevant); (b) phase slip on junction 
1 (or 2) when is relevant and slip-anti-slip pairs proliferate. 



transitions occur whose locations and properties are not 
given correctly by either the weak or strong coupHng ap- 
proaches. 

A comparison of Fig. 5 and Fig. 9 reveals that there 
is a difference between weak and strong coupling phase 
diagrams for r < Rq. In particular, the inferred phase 
boundaries between the FSC and NOR phases differ in 
these two limits. This transition is special in that both 
junctions go from superconducting to normal, but the 
transition is driven by the dynamics of just one of them 
(note that there was no direct FSC to NOR phase bound- 
ary in the naive weak-coupling phase diagram in Fig. 4). 
In the weak coupling limit, when we analyzed the su- 



perconductor to normal transition of junction one, our 
underlying assumption was that junction two was nor- 
mal. By contrast, for the same transition in the strong 
coupling case, junction two was assumed to be effectively 
superconducting. This distinction between the approxi- 
mate descriptions accounts for the difference in inferred 
phase diagrams. What is the actual behavior in this 
regime? Does it, in contrast to the other regimes, de- 
pend on the magnitudes of the Josephson couplings as 
well as the resistances? 

In Fig. 12 we indicate parts of the phase diagram for 
which weak and strong coupling analyses suggest differ- 
ent natures of the ground state. These regimes of the 
resistances would be fully superconducting (FSC) in the 
strong couphng approximation and normal (NOR) in the 
weak couphng approximation: the FSC fixed manifold 
is stable to small fugacities of the phase slips, and the 
NOR fixed manifold is stable to small Josephson cou- 
plings. This suggests that in such regimes, there should 
be a transition from NOR to FSC as the Js arc var- 
ied at a finite non-zero value of the Josephson couplings. 
Specifically, if an appropriate combination of the Joseph- 
son couplings is greater than some (resistance dependent) 
critical value then the system will be in the FSC state, 
while if this combination is less than the critical value, 
the system will be in the normal state. As such a transi- 
tion is presumably controlled by an intermediate coupling 
fixed point, it will have very different character than the 
other transitions; from now on we will refer to regimes in 
which such critical fixed points occur as simply interme- 
diate regimes. 

It is useful to remember that the original microscopic 
model had J+ = 0, so for fixed resistances in the interme- 
diate regime, on the Ji, J2 plane there will be a manifold 
below which the system fiows to the normal fixed-point, 
and above which it flows to the FSC fixed-point; this 
is the critical manifold of the FSC to NOR transition. 
Alternatively, the microscopic model could be defined in 
terms of the phase slip fugacities, Ci, C2 with C_ = 0. For 
fixed resistances in the intermediate regime, the critical 
manifold would show up here too, separating the FSC 
and the normal phase in, for fixed resistances, the Ci , C2 
plane. 

In general, an analysis of the critical behavior in the 

intermediate regime is beyond the methods of this paper, 
but we can make use of the weak and strong coupling 
limits to analyze parts of this regime: specifically, when 
the critical values of either the Josephson couplings or 
the QPS fugacities, respectively, are small. 



A. Weak Coupling Limit 

We first study the weak coupling limit. In order to 
find the critical values of Ji, J2, J+ in the intermediate 
regime, we need to analyze the effects of the non- linear 
terms in the RG flow equations (30), and, if there is in- 
deed a perturbatively accessible critical fixed point, find 
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it and the corresponding critical manifold. Truncating at 
second order, we indeed find a fixed point: 



{Jtf 

{J2f 



{R2+r- 


-RQ)iRi+R2- 


-Rq) 




rRi 




{Ri+r- 


-Rq)(Ri+R2- 


-Rq) 




rR2 




. {R2+r 


-RQ)iRi+r- 


-rQ) 


R2R1 ' 



(40) 



with an overall cutoff-dependent proportionality coeffi- 
cient having been set equal to unity when the RG equa- 
tions were first derived. As we see below, this fixed point 
can be shown to be critical provided each of the three re- 
sistance combinations in parentheses are positive. These 
factors, which we will call 



R2+r-RQ 



Ri+r-Rg 



W ': 



R1+R2 — RQ 

(41) 

are the negatives of the eigenvalues of the three couplings 
along the normal fixed manifold so that the normal phase 
is stable to small Js in this regime as indicated by the 
weak coupling phase diagram. Naively, one might have 
expected the non-linear perturbative analysis to be valid 
only when all three of these eigenvalues are small, but 
we see that in fact all that is needed is two of the three 
eigenvalues small and negative with the third being ar- 
bitrarily negative. Correspondingly, we require that all 
three of u, v, w are positive with two of them being 
small. 

By rescaling the Js appropriately, the RG flows can be 
put in a simple symmetric form in terms of u, v, w, and 
the fixed point values written as 



uw 



J* _ 

Jo — 



vw 

Vr^ 



7* — 



uv 



(42) 



rRi ^ y rR2 ^ \ R1R2 
The linearized flows around this intermediate coupling 



fixed point yield the eigenvalues which are given by 

Xi^Ai{u + v + w) (43) 

with the {Aj} being the three roots of 

A^ -I- A^ = m (44) 

in terms of the dimensionless combination of the resis- 
tances 



Auvw 



m ; 



{u + v + w)^ 



(45) 



We see immediately that for m positive, as it must 
be, there is always a unique positive eigenvalue, A+, 
which controls the growth of deviations from the criti- 
cal manifold; the two others have negative real parts and 
are hence irrelevant at the intermediate coupling critical 
fixed point. Note that if only two of u, v, w are small, 
with, say w being much larger than the other two, then 

A+ w « 1 so that A+ w 2^/uv. If all three are 

small and comparable, A+ will be of the same order but 
depend in a somewhat complicated way on their ratios. 

B. Strong Coupling Limit 

It is clear by examining the limits of validity of the 

weak coiipling expansion above that we cannot extract 
the critical behavior throughout the intermediate regime 
from this analysis. Fortunately, we can access another 
part of this regime from the strong coupling direction. 

Using the second order RG flows in terms of the fugac- 
ities of phase slips, we find a critical fixed point at 



(cr)' = 1 


( 1 R2+Ri \ ( 1 
[^Rq Y )\Rq 
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( 1 R2+Ri \ ( 1 
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Y J 
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(46) 



with Y = r(i?i + R2) + i?ii?2- As for weak coupling, manifold, defining 

it is convenient to work in terms of the negatives of the „ , 

eigenvalues of the three fugacities about the FSC fixed ^ = ~ 



v=^-Rq (47) 

W = - Rq 

with the condition for the validity of the expansion being 
that all these must be positive with at least two of them 
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small. The expansion is carried out in exactly the same 
manner as for the weak coupling limit, and the eigenval- 
ues about the intermediate coupling critical fixed point 
determined by exactly the same conditions as in Eqs. (43- 

45) , with simply u, v, w replaced by their (overbared) 
strong coupling equivalents. 

C. Superconducting-normal critical manifold 

Prom the above discussion we see that direct tran- 
sitions between the FSC and NOR phases will always 
be controlled by intermediate coupling fixed points. Al- 
though we thus cannot find the full phase boundary ex- 
actly in the intermediate region of the resistance space, 
we can use the weak and strong coupling analysis to find 
it in some regimes of the intermediate region. Eqs. (40, 

46) apply in the weak and strong Josephson coupling lim- 
its, respectively so that we can locate the phase bound- 
aries accurately in the intermediate region from the flow 
equations provided that both the bare and the fixed point 
values of the Josephson couplings are either all large or 
all small. In particular, we have found that the J* go 
to zero along certain lines in the r, R2 space which 
intersect the constant r surfaces shown in Fig. 12 at the 
points Aq, Ai, A2: our weak coupling analysis is con- 
trolled in their vicinity providing the bare Js are small. 
Analogously, the fixed point values Q vanish at points 
Bq, B\, B2 of the constant r surfaces as shown in Fig. 
12 and the strong coupling analysis is controlled in their 
vicinity provided the bare Js are large . 

The finite values of the J*'s at the fixed point on the 
critical lines has interesting implications for the phase 
boundaries in the full Rs and Js parameter space as 
sketched in Fig. 13. If we cross from FSC to NOR phase 
by changing resistances and keeping Js fixed, the exact 
location of the transition will generally depend on the val- 
ues of the Js. However, there is a whole range of small 
Js (which we can schematically denote as < J < J*) 
for which, in the second order KG approximation, this 
transition occurs exactly at the FSC to intermediate re- 
gion boundary; if we consider higher order terms in the 
RG, the location of the transition in this range will be 
modified slightly. Analogously there is a range of large 
J's for which the FSC to NOR transition happens very 
close to the intermediate region to NOR line (in strong 
coupling this occurs for < C < C*)- 

For illustrative purposes we calculate explicitly the 
phase boundary as a function of weak Ji^2 in the part of 
the intermediate regime of resistances in which the fixed 
point is at small but non-zero coupling. In particular, we 
consider the FSC to NOR transition for: 

r < 0.5, R2 = l-r + u, Ri = 1 - r + v (48) 

with u and v small, and for convenience, we set Rq = 1 
for this section. The third parameter, 

w = Ri+R2-l = l-2r + u + vf^l-2r (49) 



is generally not small. It is convenient to define rescaled 
couplings by 

which have fixed point values « y/u and K2 « y/v. 
Prom the RG flow equations, it can be seen that J+ 
rapidly approaches its nuUclinc value, ,r_^{Ki, K2), and 
then evolves slowly with the other variables. Substitut- 
ing J" for J+ in the flow equations for ii'1,2, we can flnd 
the invariant manifold on which the critical fixed point 
lies. This is parametrized by 

Kf -u[l + HKf/u)] w i^l - u[l + ln(i^:| /v)] (51) 

which has two branches of solutions; the branch with one 
of Ki or K2 larger than its fixed point value and the other 
smaller is the desired critical manifold. Note that as Ji 
increases above its fixed point value, the critical value of 
J2 decreases exponentially, and visa versa. Although we 
have taken the bare J+ = 0, even a J+ of order the fixed 
point values of the other Js will not appreciably change 
their critical values in this regime with w ^ u,v. 

Similar analysis can be done with either of the other 
pairs, u, w or v, w both small and the third of order 
unity. In these cases, however, the smallncss of the bare 
J+ means that the early stages of the renormalization will 
give rise to a non-zero value of J+ at intermediate scales 
whose value is needed to estimate the critical condition 
that relates the other Js. 

Symmetric case. Although unrealistic for the physical 
model of two junctions, it is instructive to consider the 
case in which there is a symmetry between the three su- 
perconducting components and the Josephson couplings 
linking them. In this case we take 

r = Ri=R2 = R and J+ ^ ^ = J (52) 

and the RG fiow equations become simply 

^ « J(l-2i?) + i?j2 (53) 

with 

w = u = v = 2R-l (54) 

so that the weak coupling part of the intermediate region 
occurs for R slightly bigger than ^Rq. The critical value 
of J is then simply 

Je« J*«4(i?-1) (55) 

and the RG eigenvalue controlling flows away from this 
is 

A«2(i?-i). (56) 
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FIG. 12: Intermediate coupling fixed point regions of the phase diagram. The shaded regions surrounded by a bold lino lie 
in a superconducting phase for strong-coupling and in an insulating phase for weak-coupling. For r > Rq there are no such 
regions. For intermediate J the phase boundaries will be in the shaded regions. The points Aq, Ai, A2 mark where the critical 
fixed point, J* goes to zero, and the points Bo, Bi, B2 mark where C* goes to zero corresponding to J* — > 00. Near these 
multicritical points the RG analyses in the text becomes exact. 



In the strong coupling limit, we can similarly use a single 
QPS fugacity C and write 

so that the intermediate region occurs for 

2<i?<3. (58) 

Near the upper end of this range, R shghtly less than 
|, the critical value of the fugacity is small, and the RG 
eigenvalue for deviations from criticality becomes 

X^Y--R). (59) 

Comparing the two limiting expressions for A, we see 
that, for the symmetric case, it is unlikely to get above 
a small value of order 0.2 anjrwhere in the intermediate 
region. 



VI. SYMMETRIES OF THE TWO-JUNCTION 
SYSTEM 



From the microscopic model of Fig. 2, the only obvi- 
ous symmetry — more properly a simple duality — is the 
exchange of the two junctions, R\ i?2 and J\ J2. 
Analysis presented in this section uncovers additional 
symmetries in the phase diagram of the system at zero 
temperature; indeed, in the analysis of the previous sec- 
tion we have already seen evidence of these. Here we 
will show more generally that the junction interchange 
is only one part of a larger permutation symmetry, or 
triality, that involves the interchange of all resistors r, 
R\, and R2 and the corresponding Joseplison couplings. 
We also show how the familiar weak to strong coupling 
duality of a single shunted Josephson junction"^^'^^ can 
be generalized to the two-junction system. These sym- 
metries allow one to relate in a non-trivial way many of 
the phase boundaries shown in Fig. 12. 
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1-r 



FIG. 13: Example of phase diagram in vicinity of a transition between the FSC and NOR phases for fixed r < 1/2. The critical 
manifold in the intermediate regime depends on the Josephson-coupling strengths, (a) Small J phase diagram in Ri, R2 plane 
showing FSC - NOR J = Q transition line (bold) at 7?2 = 1 — r. (b) Schematic cross-section of phase diagram along line with 
R2 = 1 ~r showing the jump in Jc suggested by the truncated second order RG analysis for crossing the phase boundary from 
ii2 < 1 — r to J?2 > 1 — The arrows indicate the RG flow of the Josephson couplings. Higher order terms in RG flows are 
likely to drive the critical Jc to zero on the line _R2 = 1 — »" (c) Three-dimensional view of the phase diagram, focusing on the 
FSC - NOR transition. The solid lines in the x-y plane mark the phase boundary between the mixed phases and the insulating 
and the FSC phases. These phase boundaries are independent of J. 



A. Permutation triality 



The two-junction system exhibits three normal phases 
and two superconducting ones. The simplest insulat- 
ing phase involves proliferations of all three kinds of 
phase slips. Conversely, the simplest superconducting 
phase, the fully-superconducting one, FSC, has none of 
the phase slips proliferating. Of the three remaining 
phases, two are normal as far as inter-lead properties 
are concerned, because of phase slips that proliferate in 
one of the two junctions. The last phase is the SC* 
phase, which is superconducting because it exhibits dis- 
sipationless lead to lead transport due to Cooper pair 
cotunneling processes. This phase, however, also has sig- 
natures of normal phases, in particular localized charges 
on the middle grain and the proliferation of QPS - anti- 



QPS pairs that decouple the phase of this grain from the 
linked superconductivity of the two leads. 

An alternative way to group the five phases is thus 
as one purely normal phase; one purely superconducting 
phase; and three mixed phases, in which part of the sys- 
tem is normal and part is superconducting. In terms of 
phase-slip fugacities, these correspond respectively to one 
phase in which all fugacities grow under the RG transfor- 
mation, one phase in which all fugacities renormalize to 
zero, and three phases in which only one of the fugacities, 
Cij C2, or grows under RG, while the remaining two 
renormalize to zero. Such grouping is very suggestive of a 
permutation symmetry of the full phase diagram in which 
the phases N1-S2, S1-N2, and SC* are transformed into 
each other, and phases FSC and NOR are invariant. In 
this section we show that such triality is indeed present 
in the low energy properties of the microscopic models 
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FIG. 14: RG flows in the intermediate region with Ri = 
0.51, R2 = 0.5001, r = 0.5. Note the typical flow pattern 

in the vicinity of the unstable critical fixed point marked by 
an asterix. (a) Projection of the RG flow trajectories on the 
Ji, J2 plane, (b) A 3D flow diagram for near-critical trajec- 
tories. 



(26) and (35) describing the system. Note that other 
systems possessing triality have been discussed earlier by 
R. Shankar in Ref.^'^; as in our case, these are non-trivial 
in some representations but easy to see in others. 

To demonstrate the triality in the original quantum 
action we consider the strong coupling representation 
of equation (35), although equivalent arguments can be 
made for the weak coupling represcntaticm described of 
equation (26). Let us begin with the mathematical for- 
mulation of this symmetry. 

The action in (35) reads: 



Z = jD[ei]jD[e2] exp -0j:\^n\^u,^R0^.+ 

dT (Cl COs(^i) + C2 COS(02) + C- cos {01 - 02)) 



where the resistance matrix is 

' r + Ri —r 
—r r + R2 



R 



(61) 



and the vector has components 0\^2- An interchange 

of the two junctions, Ri ^ R2 and (i ^ (2, will leave 
the phase diagram invariant, exchanging the two mixed 
states in which one junction is superconducting and the 
other is normal (N1-S2 and S1-N2). In (61) this in- 
terchange of junctions corresponds to transforming the 
fields, 01 02 or: 



(62) 



In terms of the new variables 0' = S ^0: 

Z = JD[0'^] J D%] exp (-a E \wn\0Z„R'0'^^ + 

dT (Cl COs(^i) + C2 COS{0'^) + C- cos {0[ - 0'2))) 

( 





where 



R' = S^RS 



r + R2 -r \ 
-r r + Ri J 

Cl = C2 C2 = Cl C- = C- 



(63) 



(64) 



This new action (63), (64), has Ri <-> R2 and Ci C2 but 
otherwise exactly the same physics with simply relabeling 
the fields 0i. 

A less trivial symmetry involves the transformation 



(65) 





leading to the action (63) with 

R' = 

Ci = 




(66) 



C'-=C2 



(60) 



This new symmetry is surprising as it swaps R2 with r. 
One way of understanding this is as a change of basis 
for the quantum phase slips. Earlier we took QPS on 
junctions one and two as a basis (schematically, we can 
label them as (1,0) and (0,1)) and considered a QPS 
dipole as their composite: (1, —1) = (1, 0) -I- (0,-1). An 
equivalent basis set, however, can be obtained by taking 
one of the QPS and the dipole as the basic objects, 
and viewing the other QPS as their composite e.g., 
(0, 1) = (1,0) -|- (—1, 1). The corresponding transforma- 
tion (65) maps phases S1-N2 and SC* into each other, 
while leaving the other ones intact. 
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Using transformations (62) and (65) one can construct 
transformations that permute any of the three resistances 
and connect any of the phases N1-S2, S1-N2, and SC* . 
The physical basis of this symmetry foUows from the 
observation that the circuits corresponding to the three 
kinds of phasc-shps arc similar; one resistor is connected 
in series to the two other resistors, which are connected 
in parallel. (The strong coupling representation we are 
using here implies starting from the FSC phase as in Sec- 
tion IVc). From the circuit diagrams in Fig. 10 we see 
the origin of the permutation symmetry: circuits associ- 
ated with all three kinds of phase slips differ only in the 
exchange of resistors. The strong coupling permutation 
triality thus generally corresponds to 



with 



1,2,- 



(67) 



where we have paired the fugacities with the correspond- 
ing resistance, so that i?+ = r, and tt is a permutation 
of the three indices. 

In the weak coupling regime, the nature of the trial- 
ity is the same: the circuits corresponding to the three 
Cooper pair tunneling events are similar with two resis- 
tors in series and a third taken out of the circuit. (Use of 
the weak coupling representation implies starting from 
the NOR phase, see Sec. 111(c).) If we now pair the 
Josephson couplings with the corresponding missing re- 
sistor in the equivalent circuits, ri = R2, r2 = Ri, and 
r_|_ = r, the permutation symmetry in the weak-coupling 
limit becomes 



T' — J 




FIG. 15: Y ^ A transformation. "Y" resistor network on left 
is mapped to A network on right via ZiZa = Z2Zb = Z^Zc = 
Z1Z2 + Z2ZS + Z3Z1. The inverse transformation (A Y) 

is ZbZc/Zl — ZaZajZi ~ ZaZb/Zz — Za + Zb + Zc- 



Cooper pair tunneling event (or a QPS) in one of the 
junctions depends on the state of the other junction (see 
Sees. III(c) and IV(c)). 

The duality transformation maps Cooper pairs into 
QPS and superconducting phases into normal ones. 
Hence, when we discuss the duality between Cooper tun- 
neling events and QPS on any given junction, we need the 
duality transformation to change the state of the other 
junction. For example, consider a Cooper pair tunnel- 
ing through junction one with junction two normal. The 
dual of this will be a QPS on junction one, with junc- 
tion two superconducting. Comparison of the effective 
shunting resistances in the two cases immediately gives 
the duality relation 



with 



1,2, 



(68) 



Ri 



Ri 



rn,2 



(69) 



where tt is again a permutation. 



B. Weak to strong coupling duality 

The similar form of the strong-coupling and weak- 
coupling representations of the quantum actions (26) and 
(35) suggests that there is a duality between the two 
regimes. 

The duality we find is a generalization of that of a 
single resistively shunted Josephson junction (see e.g.^^). 
For the single junction the duality is equivalent to the 
observation that quantum phase slips in a junction with 
shunt resistance R behave similarly, as far as their 
quantum statistical-mechanics, to Cooper pair tunneling 
events in a junction with shunt resistance R — Rq/ R. In 
the two junction problem discussed in this paper we ex- 
pect that Cooper pair tunneling events across any of the 
junctions in weak coupling should be dual to quantum 
phase slips on the same junction in strong coupling; and 
Cooper pair cotunncling processes across the two junc- 
tions should be dual to QPS dipoles on the two junctions. 
But a complication is that the effective resistance for a 



Analogous arguments give 



R2 



R2+: 



(70) 



H1H2 
R1+R2 



An alternative way of seeing the duality is to take Oi — > 
A'l and 62 ^^'2 i^i action (35). The cosine terms of 
the resulting action in terms of {A[, Aj) and those of the 
weak coupling action (26) then have the same form. If 
we compare the quadratic terms in these actions, we find 
the same duality relations (69) and (70). 

We can solve the duality relations (69) and (70) for 
f, ^1, R2: 



D _ E?2 -Rl 

R2 - Rq — 



(71) 



with Y = rRi + rR2 + RiRi- This mapping of the 
resistors to dual resistors may seem rather unintuitive, 
however Eqs. (71) coincides with the well known "Y-A" 
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FIG. 16: (a) Original circuit of Fig. 2 in weak coupling limit 
showing a "Y" resistor network . (b) Strong coupling dual of 
the circuit showing a "A" shaped network. The network in 
(b) captures the duality Eqs. (71). 



transformation of resistor networks. The Y-A transfor- 
mation is depicted in Fig. 15. By comparing the Y-A 
transformation equations in Fig. 15 we see that the du- 
ality transforms the system in Fig. 16(a) to the system 
in Fig. 16(b). In Fig. 16(a) the resistors R2 and r 
are connected in a "Y" pattern; the transformed system 
has the resistances R'q/R2, Rq/Ri and Rg/r connected 
a A pattern. 

This statement of the duahty is simple; pair-tunneling 
events (current sources) with a Y resistance network and 
resistors r, R2 (Fig. 16(a)) are dual to quantum 
phase-slips (voltage sources) with a A network of resis- 
tances i?yr, Rq/R2 (Fig. 16(b)). This is a 
simple generalization of the single junction duality. From 
Fig. 16 we see that as r — > the duality reduces to: 



C. Phase boundaries controlled by weak or strong 
coupling 

The weak-to-strong coupling duality yields a mapping 

between various of the phase boundaries in Fig 12. The 
nature of this mapping is such that weak coupling tran- 
sitions will be mapped to strong coupling ones, e.g. the 
NOR to N1-S2 boundary gets mapped into the FSC to 
S1-N2 boundary. Here NOR to N1-S2 corresponds to a 
weak coupling transition, since it involves ordering of A2 
with Ai remaining disordered on both sides of the tran- 
sition, i.e., J2 becomes relevant, while Ji and J+ stay 
irrelevant. By contrast S1-N2 to FSC is really a strong 
coupling transition because it involves J2 becoming rel- 
evant with Ji already relevant. This latter transition is 
simple in terms of the QPS fugacities, corresponding to 
^2 becoming relevant about the FSC manifold with ^1 
and C- irrelevant on both sides of the phase boundary. 

First, we map the phase boundary R2 + r = Rq via 
(71). After substituting r = Rq — R2 this yields: 



r = 



which is simply the duality of a single junction applied 

to the two uncoupled junctions, as should be expected in 
this limit in which the middle grain is macroscopic. 



R1+R2+ 



Rq-R-2 



7? ^0—^2 
^ Ri+R2-Ri/RQ 



Ri 



RQ+R1-R2+ 



R2 



Ri+R2-R:i/RQ (72) 
R2 = ^ - °- «^ 



Rq I "^'"Q-"^ ~ Ri+R2-Ri/RQ 



These apparently complicated expressions are simply the 
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boundary of the FSC phase, since: 



i+- 



R1+R2+ 



Hi -"2 



-R2 + - 



(73) 



1 R1+R2-R2/RQ _ j_ 
Rq Ri+(Rq-R2)R2/Rq Rq ' 



as shown in Fig. 17(a). As a second example, consider 
the critical Une Ri + ^^'^ = 1, which separates the FSC 
phase from the mixed phase in which junction one is nor- 
mal. The duality equations yield: 



DID, «iH2%i+H2-Hq) ~ -^Q 
^1+^2 + 



Rq—Ri 
R2 



Ri = 



R\ 



Q 



fl2(«q-Hl) , Hlfi2(fiQ-fil) 



R-L + 



= Rq 



R^+R2-Rq "^(R2 + -Ri-Rq)-R2 



^1 +^2~^Q 

R2 



R2 



R. 



Q 



, R2(Rq-Ri) , «i(fiQ--Rl) 



_ n -ffi(fli + fl2-iio) 



so that 



(74) 



(75) 



which is the condition for the phase boundary between 
the normal phase and the mixed phase in which junction 
one is superconducting, as in Fig. 17(b). 



D. Duality in the intermediate region 

In the intermediate regime of the resistance parameter 
space, the behavior under duality is more complicated. 
Since the controlling critical fixed point that determines 
the fully normal to fully superconducting phase boundary 
is at non-zero Josephson coupling in this regime, the early 
stages of the renormalization will affect the location of 
the critical manifold in the full parameter space. Thus 
duality cannot be used to locate the phase boundaries. 
Nevertheless, duality is still useful in this intermediate 
region. 

The low energy properties of the system will be given 
by the effective actions that do exhibit duality. Thus uni- 
versal properties near the transitions at pairs of points in 
resistance space should be dual even when the location 
of the transitions as functions of the Josephson couplings 
are not. In particular, as we have seen in the explicit 
perturbative calculations of the critical behavior in the 
intermediate region in the regimes in which the critical 
fixed point is at either very strong or very weak coupling, 
the critical exponents, such as the RG eigenvalue A that 
controls deviations from criticality, will be universal func- 
tions of the resistances with values on twelve-member sets 
of points being the same by the duality and the three fold 
permutation symmetry. 

For the highly symmetric case, Ri = R2 = r = R, the 
duality is simply 



R = 



R: 



Q 
3i? 



(76) 




FIG. 18: Detection of the FSC to SC* phase transition. The 
transition between FSC and SC* will induce a jump in the 
effective resistance between the leads and the grain. This can 

be observed by measuring the resistance between the lead 1 
and the normal part of the grain. The resistance measured 

by n will increase from (^"^ + "^ + 7) the FSC phase 

to in the SC* phase. A similar discontinuity in the 

resistance will also occur at other phase boundaries; see the 
discussion in Sec. VII A. 



so that there is a self-dual point at i? = at which we 
expect the eigenvalue A to attain its maximum and the 
associated correlation time exponent that controls the 
scaling of the temperature at which crossover will occur 
from critical to non-critical to be minimum 

More generally, the fact that the duality of (71) in- 
volves the combination F in a simple way, enables us to 
immediately find a self dual condition: 



Y = rRi+ rR2 + R1R2 



(77) 



When this condition is satisfied, the system will be on 
the self-dual surface. In the intermediate region, we thus 
expect the exponent A to be maximal on this surface and 
decrease in both directions away from it. On this surface, 
it will presumably vary. 



VII. DISCUSSION 

A. Relation to experiments 

We now consider the consequences of the results ob- 
tained in this paper for the two junction system shown 
in Fig. 2. 

Existence of the SC* phase. One new prediction is the 
SC* phase that is superconducting for lead to lead trans- 
port but has localized Cooper pairs on the middle grain. 
A similar phase has been discussed previously in the con- 
text of one dimensional Josephson junction arrays^^'*'^. 

To observe the difference between the SC* and the 
fully superconducting phase in the transport between 
the two leads (labeled by (f)i and (1)2 in Fig. (18)), one 
must consider the non-linear behavior, as in both phases 
there is no inter-lead resistance at zero current. But the 
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transition SC* to FSC will be characterized by a discon- 
tinuous jump in the exponent of the non-linear current- 
voltage characteristics, reflecting a change in the nature 
of the quantum phase slips in the two phases. In the SC* 
phase the system behaves essentially as one junction, and 
current is carried by lead to lead Cooper pair cotunnel- 
ing processes that are shunted by the effective resistance 
Ri+ R2- At T = 0, for small currents we thus expect 

VocT' , (78) 

where 

ai = 2{Rq/{Ri + R2) - 1) . (79) 

This form will also obtain at low temperatures and fixed 
current as long as ksT < hi /e. But at low currents for 
positive temperature we expect 

V (X T"' (80) 

(see^^). 

In the FSC phase both junctions are superconducting 
and quantum phase slips can appear in each of the junc- 
tions. The shunting resistances for QPS in junctions one 
and two are i?i + rR2/{r + R2) and R2 + rRi/{r + Ri) 
respectively, so we expect at T = and small currents 

V (x[max{T,I/e)]'^^ (81) 

with 

a2 = 2{RQ/{R'Ji, - 1) (82) 

in terms of 

KiL = max{Ri + rR2/{r + R2), R2 + rRi/{r + Ri)) . 

(83) 

Another way to distinguish the FSC and SC* phases 
is to measure resistances directly between the leads and 
the grain, as shown in Fig. 18. The effective resistances 
between the grain and the leads should jump at the tran- 
sition between the FSC and SC* phases. In order to 
measure this jump, consider adding to the circuit an ohm- 
mcter, O, measuring the resistance between the (normal) 
grain and lead one. The transition between FSC and SC* 
will be characterized by the measured resistance increas- 

...... ,„,e 

change if r is small. This occurs because in the SC* 
phase the superconductivity on the grain is effectively 
decoupled to that current cannot flow through r. 

The ohm-meter could also probe other phase transi- 
tions. For instance, in the NOR and N1-S2 phases, the 
measured resistance would be while in the S1-N2 
phase, it would be 7^7^ • 

Observation of T* . Another result of our analysis is 
the existence of a new temperature scale T* set by the 
grain level-spacing like parameter S. At high tempera- 
tures, T >> T*. the Josephson junctions are effectively 
decoupled with the dissipation set by individual shunt 



resistances Ri and R2 (sec discussion below Eq. (25)). 
At temperatures below T*, in contrast, we have a system 
of strongly coupled Josephson junctions with the dissi- 
pation determined by the whole circuit. For example, 
in the case r > Rq, the effective dissipation is the total 
shunting resistance Ri + i?2. One possible way to ob- 
serve the crossover at T* is to choose parameters so that 
r > Rq, R1.2 < Rq but Ri + R2 > Rq. For T > T* 
dissipation is then strong enough to stabilize supercon- 
ductivity on the individual junctions and we expect that 
the measured resistance of the system will decrease with 
decreasing temperature. But below T* the dissipation is 
no longer sufficient to stabilize phase coherence between 
the leads as (i?i -|- R2)/Rq > 1. At this point the phase 
slip fugacities become relevant, and we expect an upturn 
in the linear resistance as the temperature is lowered fur- 
ther. The basic reason for this is that at lower tem- 
peratures, the superconductivity is determined by longer 
length-scale fluctuations that involve less dissipation; the 
superconductivity is more vulnerable to these than the 
higher temperature more dissipative fluctuations. 

Universal vs. non-universal behavior of the resistance 
at the transition. An interesting feature of the zero tem- 
perature phase diagram, which contrasts with that of a 
single junction, is the occurrence of some of the normal to 
superconductor transitions at non- universal values of the 
total resistance . Other transitions will occur at universal 
values of the appropriate resistance. 

1. In the mixed phase S1-N2, the linear resistance of the 
whole circuit is R2 + rR\/{r + Ri) (junction one is 
superconducting, and junction two is insulating, see 
Fig. 10(b)). When this resistance becomes equal to 
the quantum of resistance Rq there is a transition into 
the superconducting state FSC. That this transition 
occurs at a universal value of the total resistance is 
not surprising: it is due to the ordering of the "last" 
non-superconducting junction in the otherwise super- 
conducting circuit. 

2. In the fully normal phase, the system has resistance 
R\ + R2. At the transition point into the super- 
conducting SC* phase Ri + R2 = Rq, so we again 
have a universal total resistance. This transition into 
the SC* phase is like a global or "long wavelength" 
one: it involves superconducting fluctuations of the 
longest lengthscale available: lead to lead cotunneling 
of Cooper pairs . 

3. At the direct transition from NOR to FSC, Ri + R2 
does not assume a universal value. For example, in the 
limit of small r the transition takes place when both 
resistances are close to Rq (see Figs. 5 and 9), so the 
total resistance will be around 2Rq at the transition. 
When r — > the two junctions are decoupled even at 
zero temperatures (sec Eq. (28)). This limit is an ex- 
ample of a "local" superconductor to normal transition 
in which the resistance per junction is equal to Rq at 
the transition point. This is the limit that has been 
extensively considered in the literature. ^'''^^ 
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Tuning the superconductor to normal transition by 
changing the Josephson couplings. Wc have shown that 
the superconductor to normal transition in a two junc- 
tion system may be tuned by changing the Josephson 
couphngs, Ji or J2, as well as by changing the shunting 
resistances i?i,2- The former may be easier to control in 
experiments as demonstrated recently in Refs.^. 

Non- Universality of the critical exponents. In Section 
VI wc showed that the transition between the fully su- 
perconducting and fully normal phases is controlled by 
a fixed point at intermediate values of the Josephson 
couplings. The critical exponents of this transition arc 
non-universal and vary continuously as the three resis- 
tances in the system change. Non-universality of the 
critical exponents at superconductor-normal transitions 

in the presence of dissipation has also been discussed in 
Refs.20.21,56 

Symmetries of the two junction system. In Section 
VI we discussed the rich symmetries of the two-junction 
system. In addition to the usual weak-strong coupling 
duality'^^'^^ it exhibits a permutation-triality. Exchang- 
ing the three resistances Ri, R2, r leaves the action and 
the phase diagram essentially unchanged. These sym- 
metries provide a powerful tool for studying the two- 
junction system; one need only investigate one corner 
of the phase diagram to be able to construct it in its en- 
tirety. The boundaries of the region in which there is an 
intermediate coupling fixed point (see Fig. 17), can be 
found from the triality and weak-strong duality transfor- 
mations. 



B. Broader relevance and open questions 

The results obtained in this paper should provide hints 
that may help understand other superconductor to nor- 
mal transitions, such as in thin wires^^'^^ and in films^. 
It is often conjectured that such transitions can be de- 
scribed in terms of models of resistively shunted Joseph- 
son junctions. For example, in wires one might perhaps 
think of segments of wire of length ^0 (i-e. the super- 
conducting coherence length or phase slip core size) as 
individual grains. Then to estimate the crossover tem- 
perature analogous to our T* one could take both R and 
r of order the normal state resistance of a single segment. 
This would yield a superfluid-to-normal relaxation rate 
that is of the order of Tc. The crossover temperature, 
T* , is related to the energy level separation parameter 
S in such a segment of wire of length ^o- Using dirty 
limit exressions = 1.87lD/Co, R = i,^/{e^N^DA) and 
5 = (A^o^Co)^^j in terms of D, the diffusion coefficient, 
A^o, the density of states per unit volume, and A, the 
wire's cross-section, we find T* « T^. So at all tem- 
peratures one should consider the effects of interactions 
between the effective "Josephson junctions" that link the 
"grains"; i.e., effects analogous to those discussed in this 
paper. 

One possibility is that for wires much longer than 



the superconductor to normal transition will be deter- 
mined not by the resistance per coherence length, but by 
the total normal state resistance. Such behavior has been 
observed recently in experiments of Bezryadin et.al.^'^ 
where wires as long as fifteen times ^0 had a normal to 
superconductor transition when their total normal state 
resistance was close to Rq (sec however [15,54,55]). 

There is, however, another effect that must be con- 
sidered in the long-wire regime. When normal metallic 
wires are long enough that their resistance is of order 
^ = 4i?Q, localization effects start to be important at 
low temperatures, specifically below the temperature at 
which the inelastic mean-free path of the normal elec- 
trons is of order the length over which the wire has re- 
sistance of order 4Rq. It is thus not clear that there is a 
regime in which the dissipative effects discussed here can 
affect the superconductivity without localization effects 
also becoming important. At least naively, however, sec- 
tions of length ^0 can not have resistance Rq for T < T^, 
and the inelastic scattering length is smaller than the 
coherence length near Tc. Thus there may well be tem- 
perature regimes in which these collective effects are im- 
portant but localization effects not. This clearly requires 
substantial further thought. Alternate geometries, such 
as configurations with a metal layer underlying the super- 
conducting wire, may be the best candidates for avoiding 
some of these complications. 

In the previous subsection we discussed the possibility 
of a surprising phenomenon in the two-junction system: a 
minimum of the resistance at a crossover temperature T* 
with an upturn at lower temperatures. Qualitatively sim- 
ilar behavior has already been observed in experiments 
on Josephson junction arrays and superconducting films. 
It is likely that the disorder plays an important role in 
such systems — especially in granular films such as InO.'^ 
Close to superconductor to normal transitions in disor- 
dered materials, the behavior may be dominated by weak 
links that involve connections via mesoscopic size grains. 
As the temperature is lowered below the local T*, the 
effective dissipation shunting these links will change in a 
manner analogous to that of the pair of junctions in se- 
ries through a small grain discussed in this paper. This 
could potentially account for the observed saturation of 
the resistance at low temperatures in systems that would 
appear to be becoming superconducting on the basis of 
their behavior at higher temperatures. Understanding of 
such systems would benefit from generalizing the analy- 
sis of the two junction system presented here to arrays of 
superconducting grains and Josephson junctions in both 
one and two dimensions. 

An important issue that we have not addressed is the 
microscopic nature of the charge relaxation between nor- 
mal and superconducting fluids that we have introduced 
phcnomcnologically. We have assumed that at low fre- 
quencies this is ohmic even in the limit of zero temper- 
ature, but even if this is indeed the case, r should cer- 
tainly depend on details of the experimental system. If, 
in fact, the relaxation is subohmic or superohmic in the 
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low temperature limit, this will be roughly equivalent to 
the r — > oo of r — ^ cases discussed here. However taking 
into account charge quantization effects on the super-to- 
normal fluid relaxation and the role of quasiparticles and 
their non-conservation may lead to qualitatively new ef- 
fects. One question that must be considered is whether 
there will be enough low energy excitations on scales be- 
low T* to give rise to the dissipative effects that are cru- 
cial for the logarithmic dependence of the effective action 
of quantum phase slips on temperature. We leave these 
issues for future research. 



VIII. SUMMARY 

In this paper we have analyzed Cooper pair tunneling 
between two macroscopic leads via a mesoscopic super- 
conducting grain in the presence of ohmic dissipation. 
We treated this system in terms of a two fluid descrip- 
tion to the grain by effectively splitting it into normal 
and superconducting parts with capacitative and gal- 
vanic couplings between the Cooper pairs and normal 
electrons. A phenomenological ohmic resistance, r, was 
introduced to describe the charge relaxation between the 
superconducting and normal parts of the grain. The cor- 
responding microscopic Hamiltonian was used to derive 
the quantum action in terms of which the analysis was 
carried out. We showed that there is a new tempera- 
ture scale T* that separates two very different regimes. 
For macroscopic grains, T* = 0, so that the system is 
always in the high temperature regime in which the two 
junctions are decoupled. In contrast, for small grains at 
temperatures below T* there is strong coupling between 
the junctions and the system can be described by a two 
component sine-Gordon model. We analyzed this model 
in the limit of weak Josephson coupling and showed that 
it leads to a rich quantum phase diagram with two super- 
conducting and three non-superconducting phases. The 
most surprising result is the appearance of a novel super- 
conducting phase, SC* that has loc;alized Cooper pairs 
on the grain but phase coherence between the leads due 
to Cooper pair cotunneling processes. 

The limit of strong Josephson coupling was studied 
using a dual two component sine-Gordon model. Simple 
circuit theory for the two-junction system enabled us to 
derive the phase diagram for both the weak and strong 
Josephson coupling limits. In contrast to the single- 
junction case, we demonstrated that the strong and weak 
coupling analysis predict different locations of the transi- 
tion between the fully superconducting and fully normal 
phases implying the existence of an intermediate coupling 
fixed point controlling this transition. We analyzed the 
renormalization group flows in this intermediate regime 
and found non-universal critical behavior with the expo- 
nents depending continuously on the resistances involved, 
The rich symmetries of the two component sine-Gordon 
model include weak to strong coupling duality and per- 
mutation triality of the shunting resistors i?i,2 and relax- 



ation resistance r. 

Experimental implications of our model, including the 
crossover temperature T*, the identification of the novel 
superconducting phase SC*, and the lack of universality 
of the measured resistance at the superconductor to nor- 
mal transition were discussed briefly. Finally, we noted 
that our results may be useful for understanding some 
of the puzzling properties of superconductor to normal 
transitions in thin wires and films. 
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APPENDIX A: MICROSCOPIC MODEL 

1. Microscopic model for a two-fluid network 

In this Appendix we provide the derivation of several 
important results used in Sec. II B. For generality, the 
first part of our analysis is not restricted to the system 
shown in Fig. 2, but applies to any two-fluid network. 
The network consists of superconducting islands (which 
may be electrodes or grains). Each island i in this net- 
work is assumed to have part of its charge in the form 
of superconducting Cooper pairs, Qsi, and part of the 
charge, Qni, in the form of normal fluid. The Hamilto- 
nian of the system consists of three pieces: 

H{QNi,Qsi,<Pi,i^i) = Hq + Hj + Hdis. (Al) 

The charging part Hq is given by Eq. (6), with Kij de- 
fined as in equation (5). The Josephson energy of the 
Cooper pair tunneling between the grains is 

= - 5 J^3 C0S(.^, (A2) 

Dissipation between the islands, as well as charge relax- 
ation between the Cooper pairs and normal fluid inside 
the islands, is described using the Caldeira-Leggett heat 
bath model (see discussion in Sees. II A, II B) with resis- 
tances Rij and ri respectively. 

'Hdis — \ J2ij T^bath{Rij — 2V'j) 

(A3) 

The commutation relations between charges and phases 
are given by equation (7). Note that the Heisenberg equa- 
tions of motion on and tpi correctly reproduce Joseph- 
son relations as in Eqs. (8). 
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We use the Hamiltonian (Al) and the commutation 
relations (7) to construct the imaginary time quantum 
action 

Z = J VQNiVQsiVcPiVtIJi 
exp (2ie J2i Jo dr Qsi + ie J2t Jo dr Qni ipi (A4) 
- Jl^dTn{QNi,Qsi,<l>ui'i)) I 



We remind the reader that in the presence of ohmic dissi- 
pation the phase variables, (pi and ipi, should be periodic 
at r = and t = jS (no phase twists by multiples of 27r 
are allowed). 

After integrating out Qni and Qsi in (A4) we find 



Z = J Vcji^Vipi exp {-Sq - Sj - Sdts) 

= -\Y.ijJo dTJijCOs{(t>i-(t>j) 



where the matrices M satisfy the equation 
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(A6) 



where we defined (Eq. 5) 



(A7) 



NiOij ■ 



Dsi, are the level spacings of the island i, and Cij 
is the capacitance of the island network. 

In mesoscopic grains, level spacings are already much 
smaller than the electrostatic capacitances, and this con- 
dition is even better satisfied in macroscopic electrodes. 
Hence, we can expand (A6) in Ds,n- It is useful to point 
out that this approximation does not require that every 
Ds.Ni smaller than any island of the C~^^ matrix, but 
only that Ds,Ni is smaller than C^^^ . Hence, this expan- 
sion can be applied even when we have a combination of 
macroscopic electrodes and mesoscopic grains. We ob- 
tain 



Msij 

MNSij 



Dsi + 



where 



Dsi + Dni 



Dsi + 



(A8) 



Dsi + Dni 



Vi = 



Dsi 



Dsi + Dm 

Therefore, we can use the following simple expression 



(A9) 



■i (Dsi+DNi) 



Sq = Jo dr (2(5i)7 E 



(AlO) 



The first term in (AlO) tends to equilibrate the normal 
and the superconducting fluids, by introducing an ener- 
getic penalty for having different chemical potentials. For 
macroscopic grains level spacings are zero, so this term 
requires (j) = 2ip, which is the case considered in the lit- 
erature previously. The second term in (AlO) describes 
the usual Coulomb interaction between the islands, but 
the potential on each island is now give by the weighted 
average of the potentials of the two fluids: 



Vi = {siVsi + 'HiVNi)- 



(All) 



2. Equations of motion 

As a consistency check on the quantum action (A5), it 
is useful to show that its equations of motion reproduce 
the familiar equations of electrodynamics. After taking 
functional derivatives of (A5) with respect to (pi and V'j 
and analytically continuing into real time, we have 
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(2i)2 Ej Msij '<i>j + Ej MsNiji>j - Jijsin{4>i - + ji^i^i - '2ipi) = 



(A12) 



From (5), (8), and (A6) we have 

j 3 

= ;5Z^Jivy'/'j + ^l]Af;vyV'j- (A13) 

3 3 

Eqs. (A12) may be written then as 

- — 2^ J^jSin{(j)i - H = 



dt 
dt 



VNi — VnJ _ Vsi — VNi 



= 0. (A14) 



These are the usual charge conservation equations; the 
Josephson form of the Cooper pair tunneUng current, and 
Ohm's laws for the normal currents and the "conversion 
currents" between the Cooper pairs and the normal fluid. 



Two leads Josephson coupled via a mesoscopic 
superconducting grain 



We now apply our general discussion from Sec. A 2 to 

the system shown in Fig. 1, a single mesoscopic grain 
between two superconducting electrodes. We assume 
that the electrodes are sufficiently large, so the super- 
conducting and normal fluids are perfectly coupled in 
them, 01 — 2tpi and 02 = 2V'2. From equation (AlO) 
the charging part of our system can be written as 



Sq 



2{2ey 



■ f dTY,XiC%X3, (A15) 



where = (^i, ^2, (t>g, '2'ipg), and 



J 



/ Cll Cl2 ClgSg C^lgTJg 

Cl2 C22 C2gSg C'2gr?g 

ClgSg C2gSg CggSg + Cq CggSgT^g " Cq 

\ClgVg C2gVg CggSgr]g-CQ CggT^l + CQ 



(A16) 



where C, 
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D 



Sg 



D 



Ng, 



D]yg/{Dsg + Dpfg), and rjg — Dsg/ {Dgg + DNg)- It is convenient to change variables 



to the phase differences and the center of mass phase defined as 



Ai 
A2 

A. 
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Ig I , C22 + f '12 + C2g 

+ 7^ ' 

'^tot 



Clg + Ciu + ( 



Clg + ( '20 + C, 



2a 
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Cu 



Cu 



Vg2i;g, (A17) 



with Ctot = Cii + 2Ci2 + C22 + 2Cig + 2C2g + Cgg. We 
have 

where the indices a and /? are summed over 1,2, and g. 
It is useful to observe that the center of mass phase, 
is decoupled from the phase differences in (A18) and can 
be integrated out in the partition function. 



We do not discuss the most general case of the capaci- 
tance matrix , but concentrate on the situation when 
the dominant capacitances are the mutual capacitances 
between the electrode one and the grain, C\ , and the elec- 
trode two and the grain, C2. This corresponds to taking 
Cii = Ci + ACi, C12 = 0, Clg = —Ci, C22 = C2 + AC2, 
C2g = -C2, and Cgg = Ci + C2 + ACg. After some 



29 



Superconducting 
grain 

r 



Normal Grain 



Normal Grain 



9, SC 



1/D 




Q + Q 

SC N — 



FIG. 19: The two fluid model description of a free-standing, 
mesoscopic, superconducting grain. The grain is split into two 
grains, a suporconducting-fluid grain, which contains Cooper 
pairs, and a normal-fluid grain, which contains the normal 
electrons. Normal electrons can become superconducting by 
flowing through r. The potential on the grains is given by a 
sum of the electrical potential, {Qn+Qsc)/C, and a chemical 
contribution, DnQn, and DsQsc- The finite level spacings 
are modeled as capacitors with capacitances 1/Dn, 1/Ds- 



straightforward manipulations, we get 



+C2(A2 + Vg^gf + CqAI + Ctot^^) . 



(A19) 



4. Circuit-theory approach to the two-fluid model 

We can gain more intuition about the analysis pre- 
sented in Section IIA by considering effective circuits for 
the island network. As a first example, let us take a free 
standing grain. The electrochemical potentials for the 
normal and superconducting electrons on the grain can 
be written in the form 



Vg,sc = ^^'^ + DsQsc. 



(A20) 



Here C is the capacitance of the grain relative to the 
ground, and the Dj's are the inverses of the correspond- 
ing compressibilities. Eq. (A20) describes the electrical 
system in Fig. 19. In addition to C, there are two more 
"effective" capacitors, l/Dsc, which describe the 

extra potential drop produced by the level spacings in 
each part of the grain. As can be seen in Fig. 19, the 
charge on the capacitor C has to be equal to the total 
charge on the grain, Qn + Qsc- 

The electro-chemical potentials in (A20) yield the 
charging part of the Hamiltonian: 



T^Q - WjiQ 



2C 



!SC 



h^sQsc- 



(A21) 




Superconductirtg Grain 



FIG. 20: The effective circuit of the two Josephsori junc- 
tion system. The mesoscopic grain is cormocted to the leads 
through Josephson junctions and resistors. It also interacts 
capacitatively with the leads. This interaction is modeled 
by the capacitors Ci, Ci which connect to the two parts of 
the grain through additional capacitors, 1/Dm, 1/Ds- The 
additional capacitors account for the finite level spacings in 
the grain. The "bare" electrical potential on the grain (the 
electro-chemical potential without the level-spacing contribu- 
tion) is given by Vb, as noted in this figure. 



Prom here on we could proceed along the lines of Ap- 
pendix Al to obtain the action for this circuit. 

The general principal behind Eqs. (A20) is that the 
potential on each island consists of a sum of the electrical 
contribution, Ve, due to Coulomb interactions, and the 
level spacing contribution: 



Vg^ n=-Ve + DnQn 

Vg, SC = + DscQsC- 



(A22) 



If we construct a circuit for an island network, Eqs. (A22) 
indicate that we need to put the extra effective capaci- 
tors, l/Djvi, 1/Dsi, between the point at which a macro- 
scopic island would be, and the normal and supercon- 
ducting grains, respectively. Let us demonstrate this by 
constructing the effective circuit of the two-junction sys- 
tem. 

The two junction system consists of a mesoscopic su- 
perconducting grain situated between two macroscopic 
superconducting leads (Fig. 20). The capacitors Ci and 
C2 describe the "bare" interaction between the leads and 
the grain. In addition to them, there are also the effective 
capacitors l/D^ and 1/Ds, which describe, respectively, 
the level spacing of the normal part and the supercon- 
ducting part of the mesoscopic grain (Fig. 20). These 
capacitors connect the point Vq, at which a macroscopic 
grain would have been, to the normal and superconduct- 
ing parts of the mesoscopic grain. 

The electrostatic part of the Hamiltonian of the two- 
junction system as shown in Fig. 20 is given by 



— 2Ci ^1 + 2C2 '^2 



with the constraint 



Q1+Q2 + QN + Qsc = 0. 



(A23) 



(A24) 
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This constraint merely reflects the fact that the capac- 
itors 1/Dsc and l/Dj^ are not real capacitors, but an 
electrical analogy to the effects of the level spacings 
in the mesoscopic grain. The charge on the grain is 
~Qi — Q2 (where the minus sign is due to the conven- 
tion in Fig. 20), and it is split into a superconducting 
part, Qsc, and a normal part, Qn. In turn, Qsc and 
Qn increase the electro-chemical potential on the grain, 
which is taken into account using the fictitious capaci- 
tors, l/Dsc, I/Dn- 

We can use the constraint (A24) to eliminate the 
charge of the normal grain: 



(A25) 



(Qi + Q2 + Qscf + \5scQlc- 

One can now proceed by defining the phases ^1, 02 and 
(j)g, which obey the commutation relations 

=-2ie 

[(32, M = -2ie (A26) 
[Qsc,<f>g] = -2ie. 
and following steps presented in Appendix A 1. 



APPENDIX B: LOW TEMPERATURE 
DISSIPATION 

In the discussion in Sec. II we introduced the normal 
fluid of gapless quasiparticles as the origin of the dissipa- 
tion for the junctions. This is not, however, a unique way 
of getting dissipation, including its ohmic variety. From 
the various possibilities, let us mention exciting electro- 
magnetic waves in the environment by fluctuations of the 
voltage and charge on the junctions. A well studied ex- 
ample is a junction connected to an LC line (see Fig. 
21(a)). Sudden changes of the voltage in the junction 
excites plasmons, which carry the energy off to infinity 
(away from the junction) leading to dissipation. It can 
be described using effective impedance formalism:^'' 



5*. = /3^Re 



R 



Q 



(Bl) 



where A(f> is the phase difference across the junction 
and Z{ll!) is the impedance of the environment seen 
by the junction. In the case of an infinite LC line 
Z = [Lq/Co)^/"^, where Lq and Co are inductance and ca- 
pacitance per unit length respectively, so we arrive at the 
Caldeira-Leggctt type ohmic dissipation given in equa- 
tion (14). For the system considered in this paper (see 
Fig. 1) such LC line (or its analogues) may come from 
the edges of the electrodes or the connecting wires. The 
crucial observation is that different Cooper pair tunneling 
processes (between the two electrodes and the grain, and 
the cotunneling process) should in general excite differ- 
ent electromagnetic waves. This can be seen from the 



schematic circuit shown in Fig. 21(b). The effective 
transmission lines in the figure give rise to three differ- 
ent resistors, which are related to Ra, and r from 
the model in Eq. 26, as discussed in the caption of Fig. 
21. It is useful to point out that in this model we can 
relax the assumption of the small size of the grain, since 
the electromagnetic interactions discussed here, as well 
as the Josephson couplings, are present at all tempera- 
tures. In that case wires can be connected to each of 
the superconductors separately, allowing direct measure- 
ment of the rich phase diagram discussed in the bulk of 
the text, and effects related to charge discreetness are ex- 
pected to be less significant. Dissipation due to other low 
energy degrees of freedom in the system^® is also possible. 

It is worth emphasizing that the precise form of the 
quantum model for dissipation depends crucially on its 
nature. A common choice of the Caldeira-Legget ohmic 
heat bath model comes from the fact that it is the sim- 
plest quantum model consistent with the classical equa- 
tions of motion. One expects that many effects of the 
dissipation would be at least qualitatively independent of 
its nature,^® although considerable differences may also 
be present. 



APPENDIX C: FREQUENCY SHELL RG 

In order to find the phase diagram of a resistively 
shunted Josephson junction in the weak or strong cou- 
pling regime, it is best to employ a frequency shell RG. 
Generally we start with a sine-Gordon partition function 
such as: 



exp 



duj 
2^' 



R\uj\ 



C 



dr- cos{9) 
a 



2ttRq 

(CI) 

where we have taken the T — » limit and changed the lo 
sums into integrals. It is useful to redefine the amplitude 

of the anharmonic term using the short time cut-off a ^ 



— , wher 



is the plasma frequency of the junction. 



The sharp high frequency cut-off we use is 



A; 



a 



(C2) 



As is well known^'*, the partition function (CI) is also 
the partition function of an interacting Coulomb gas in 
one dimension, with fugacity C and interaction energy: 



E, 



ij (t) 



= —2aia 



Rq 
R 



log 



(C3) 



where cr,; is the charge of the i'th particle. In the 
weak coupling limit (III A) the "particles" are pair-tunnel 
events, and C = J. In the strong coupling limit (IV) the 
"particles" are quantum phase-slips. 

Presently we would like to integrate out the " fast" de- 
grees of freedom associated with the field 6 in (CI). This 
has the physical meaning of reducing the frequency cutoff 
A, and can be thought of as increasing the effective size 
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a. 



b. 




FIG. 21: An infinite transmission lino as a source of dissipation, (a) A fine of coils and capacitors, Lo, Co, has an effective 
real impedance Z — \J Lo / Co at low frequencies. Tfiis line could describe electromagnetic modes that are excited by tunneling 
of Cooper-pairs across junction J. (b) The two-junction system may have more than one transmission line. These lines are 
reduced in the figure to the effective impedances Za, Zb, Zc- The schematic circuit shown has three effective shunt resistors, 
as in the model in Eq. (26). This configuration of effective impedances is the same as the "A" resistors network shown in 
Fig. 15. This will translate (using the Y-A transformation, Fig. 15) to the model in Eq. (26), with R\ = za+^^'+Zc ' ~ 

Zb Zg ^ _ Za Zb 
Za+Zb-\-Zc ' Za+Zb+Zc ' 



of a particle, and eliminating all particle - anti-particle 
pairs whose separation is lower than this size. We write 
the action (CI) as: 



l\uj\<A-dA 2tt 2ttRq 



which we can write as: 



duQ2 



= Z<C{l + JdT 



2a 



= C/D[e<]exp(-i; 



du o2 RW\ 
w|<A-dA 2nRQ 



From this we obtain the RG flow: 
^ R A 



(C4) 



(C5) 



We still need to restore the variables to their original 
scale so that Eq. (C2) is fulfilled. Since A ^ A - dA, 
a ^ a + da. This leads to: 



C _^ C a+da (l _ RQ_dA'\ 
a a+da a \^ RAJ 

C f 1 _ RQ_dA I dA\ 
a+da \ R A A J 



(C6) 



C(> 



Rq dA 
R A 



+ 



dl -"-dA 



CI 



dA' 
Ay 

Rq 
R 



(C7) 



Where the minus sign on the LHS of the middle equation 
denotes the fact that A is decreasing, and dl = — dlogA 
with I the differential logarithmic flow scale-parameter, 

A = Aoe-'. 

Prom (C7) we see that when R > Rq ^ is relevant 
and the particles proliferate. When R < 1 the opposite 
happens: all particles form dipoles that disappear when 
the scale increases. 



APPENDIX D: COULOMB GAS 
REPRESENTATION OF THE WEAK COUPLING 
LIMIT 

In this Sec. D 2 of this appendix we will derive the 
RG equations to second order of the two-component sine- 
Gordon model of Eq. (26). Before doing that, we will de- 
rive the Coulomb gas representation of this model (Sec. 
Dl). This representation shows that the model (26) 
describes a gas of interacting pair-tunnel events. The 
Coulomb gas representation makes it conceptually easier 
to derive the second order RG equations. In Sec. D 3 
we use the Coulomb gas description to demonstrate how 
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proliferated pair-tunnel events (or their strong coupling 
counterparts, phase slips) screen other events. This adds 
to the discussion of the mixed phases in Sec. IIIB and 
Sec. IV B. 



of J's leads to products of terms like cxp(i((/)g — 
which are translation operators for the charge-difference 
between the middle grain and lead 1. 



1. Coulomb gas representation 

To analyze the two junction system we use the map- 
ping of the partition function (26) to a partition function 
of a Coulomb gas. The starting point for this investiga- 
tion is the free energy that appears in Eq. (26): 



J 27r 2-K 



A2 



A2 M(l+- 



+2A1A2 



R1+R2 + 



+ / ^ ( Ji COS Ai + J2 cos A2 + J+ cos(Ai + A2)) 

(Dl) 

where we have redefined the anharmonic terms by a fac- 
tor of a ~ i^p^- 

The first step is to make use of the weak coupling state- 
ment, Ji, 2 ^ 1 and expand the exponent in a power law 
in the J's. Following that, an integration over the fields 
Ai^2 reduces action (Dl) to a partition function of an 
interacting gas with two kinds of charges, /x = 1, 2 cor- 
responding to exp («A^): 



exp -i/dndr2 E pS>S<:;=J 

\ /il,Al2 = l . 



(D2) 



where 



P(t) ~ S '^i^(r-Ti) is the density of the gas, and 

i—l 

ai is the charge of the i'th particle. Ji and J2 play the 
role of fugacities for the two types of gas particles. 
The interaction energies are: 



eI'I^ = -2 



(Rl+r) 



log I - I 



e: 



(22) 

(r) 



=2 



Rq 



Rq 

log I - 1 

Q => I a I 

log 1^1' 



,iR2±r 



(D3) 



where E"^^ is the interaction energy between a typc-i par- 
ticles and type-j particles. As we can see, the log diver- 
gent interactions impose the neutrality condition satisfied 
in (D2). Notice that gas particles of type-1 and type-2 
of the same charge actually attract. 

The meaning of each of the particles is very simple; a 
particle represents a cooper-pair tunneling event through 
the corresponding junction (see Fig. 3).'^'^ To see this 
recall that, for example, Ai = 0g — (/>! and the (/)'s are 
conjugate to the number of cooper-pairs on the corre- 
sponding grain or leads, hence the expansion in powers 



2. Two component gas RG 

To find the phase diagram of the two-junction sys- 
tem, we need to use both the mapping to a Coulomb gas 
from D 1, and the angular-frequency RG as described in 
Appendix C, extended to include second order contribu- 
tions. In order to make the discussion general, we will 
treat the following form of the action (26): 

Z = JVA1VA2 e-^'^-^-' 

= /o'^f (-JlCOs(Ai)- J2C0S(A2)- J+cos(Ai+A2)) 

(D4) 

with the capacitativc part omitted, the sums over a;„ ap- 
proximated by integrals, and the upper frequency cutoff 
A = tt/o introduced. The first order contributions to the 
RG flow equations come from the terms linear in the Js 
and exactly following Appendix C, but using 



(D5) 



27r 

we get the following first-order RG equations: 

^ = J, (1 - ^G-) 
^=J2(l-^G2-2^) (D6) 
^ = J+{^-^ (Gr/ + G2-2^ + 2Gr2^)) . 

The interactions between the two gas components give 
rise to second-order contributions to the RG equations. 
In addition, the second order terms in the power-law ex- 
pansion in J's of (D4) produce corrections to the plasma 
frequency of the problem and other irrelevant operators. 

First we will demonstrate how to derive all second or- 
der contributions to the flow equations by deriving one 
such contribution. Let us consider an example for a sec- 
ond order term. Then we will proceed to derive the unim- 
portant plasma-frcqucncy corrections. 

Consider the term that results from the product of a 
Ji and J2 first order terms In a power law expansion of 
D4: 

Z =... + ! D[A^]D[A^]exp (- / f^AtcA) 
/ dnj (iT2^^(e*(^i(-i)+^2(-2)) + e'(^i(-i)~^2(-2)) 
+C.C.)) + ... 

(D7) 

At this point we would like to integrate out the fast 
modes of the fields Ai, A2 as in (C4). But we need to 

be careful since if ti = T2 the suppression resulting from 
the contraction of the fast modes is not the product of 
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the two factors obtained from the first order terms in Ji 
and J2: 

(cos(Ai+A2))^>_^> = 
(-i L>A < + A>)2 >) cos (Ai + A2) 

= exp (-i ^ (Gn + + 2Gr2^)) cos (Ai + A^) 

(D8) 

The rcnormahzcd second order term, as in (D7), will 
only contain the self interaction of a Cooper-pair tun- 
neling event in junction 1 and 2, completely dropping 
the exp(2G^2^). This difference produces a J+ renor- 
malization term. To calculate this term we first need to 
separate the contribution to the partition function that 
comes from the term (D7) to same time and different 
time contributions. Define r = x = T2 — ti and 

write: 



/ cLtx J dT2 = J dr J dx = 

/ iix I >a+da d^ + U^Il^Ka^^^ + Jd^IaKl^Ka+dad^ 

(D9) 

The first integral is unaltered in the RG step (except 
for the influence on Ji of the integration of fast modes as 
in first order) and can be re-exponentiated. The second 
term represents two gas-particles of type 1 and 2 occur- 
ring at the same time, with the resolution of this RG 
step. This event too should be re-exponentiated since it 
is also obtained as a second order event in the renormal- 
izcd variables. However, as pointed out in the previous 
paragraph, there is a discrepancy in the RG suppression 
coming from the fast modes contraction. Hence we write: 



- n^JrJ^ exp (-1 t (Gri^ + G^^)) 

(cos (Ai - A2 (i + (exp (-1 /^^^, ti-^)Gi,') i)) (DfO) 
+ COS (Ai(.) + A^m) (f + (exp {-U.>A ti^)Gi2) - l))) 
= / ((1 + ^G-^) cos (A,(.) - A2(.)) + (1 - ^G-^) cos (A,(.) + A,,^,)) . 



The second term in each of the brackets multiplying the 
cos terms are the corrections that feed into J_ cos(Ai — 
A2) and J+ cos(Ai + A2). 

Eq. (DIO) leads to an additional J1J2 in + . The 
same could be done to second order terms that are prod- 
ucts of J+ and Ji^2- For instance, in the case of J_|_ the 
complete flow equation to second order would be: 

^ = ^+ (1 - 2^ (Gn + G2-2^ + 2Gr2^)) ^^^^^ 
I 



The same equation with the sign of G^2^ reversed applies 
to J_. Similarly, the flow equation for Ji (J2) would have 
a contribution proportional to J2J+ {JiJ+) and J2J+ 
(JiJ+). 

The third term in Eq. (D9) leads to terms in the ac- 
tion proportional to and hence are unimportant. To 
see this we will use two examples that exhaust all possi- 
bilities. As the first case, let us look at the term in the 
second order expansion (D7) proportional to Jf : 



Z=... + J D[A]exp(^- J ^A^GA^ J dn j (ir2^ (e*(^i(-i>+^i e^^^^f-i'-^^ c.c.) +.. . (D12) 

We ignore the first term in the brackets as it implies a very costly configuration of two particles very close to each 
other, and concentrate on the second: 

/ D[A]exp(- I — A^GA^ ( dr ( ^.t^ fe'('^i(-+-/=)-^i c.c.) (D13) 

J V i 27r J J Ja<\x\<a+da ^ ^ 



34 



Now: 



(D15) 



f^i 1^ (e'(^i(-+-/2)-Ai(x-x/2)) +C.C.) 

Jo<|x|<o+aa V / 

Rj2(^e''^iM« + e-''^iM'') (D14) 

but this can be re-exponentiated to give a correction piece for the action: 

AS= [ drda^Al,^. = [ ^da^co^Aj, 

which is just a uj^ contribution that rcnormahzcs the plasma frequency of the model. 

A more complicated case would be considering again the term that mixes the two components of the gas: 

Z = ... + jD[A]exp{-J^A'^GA) 

Here we use the following derivation: 

/a<|x|<a+da (^'^'^''"+"^'''^^''"""^''^ + ^-C-) ^ da ^e^^'^i M ^ ^e*(Ai, M-^s, m)« + e-'^^^i. M'^^^. M^") +c.c) 

« 2rfaC0S (Ai (^) + A2 (r)) (2 - ^ (Al, {r) - A2, (r) ' 



(D16) 



This results then in the introduction of a term: 



/Ji J 2 
drda—^ cos (Ai (^) + A2 (r)) 



(D17) 



(D18) 



once again, proportional to w^, and hence, unimportant. Similarly: 



Ja<\x\<a+da V V V J J 



2da 



cos (Al (^) - A2 (^)) (2-3^ (Al, (^) + A2, (^)) ' 



Giving: 



/■ , , J1J2 ^ (Ai,(^) + A2,(^)) 

A5 = y drda^ cos (Ai (,) - A2 (,)) ^ ^ ^ 



(D19) 



(D20) 



This exhausts all second order contributions to the RG flow equations. 



r 



3. Screening of pair tunneling events be understood as screening of charges of one type by pro- 

liferated charges of the other type. 

Let us start by considering the case of J2 being rele- 
When discussing the phase diagram obtained from the vant. The simplest approach to the strong coupling limit 
RG flow Eqs. (30) in Section III, we had to account is to assume that A2 = 0. This can be done because 
to parts of the phase diagram in which one of the three if J2 > 1; then the term J2COSA2 in the action (26) 

constrains A2 to be 0. The kinetic part of the action 



Joscphson couplings is relevant; then we used the pro- 
cedure of setting the respective phase- difference variable becomes: 
to zero, which is equivalent to starting from a new fixed 
point. Here we review this step and show that in the lan- 
guage of the Coulomb gas analogy, this procedure may 



Si 



-I 



^^^2 

27r 277 ^R^+R^ + Ei^ 



(D21) 
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and hence the flow for J\ (or J+) would become: 



and hence the flow for J\ (or J2) would become: 



dJi 

HT 



Ji 1 



R1 + U2 



R1R2 



(D22) 



shifting the phase boundary S1-N2 and the FSC phase 
to: 



1+32 



Q 



(D23) 



The above calculation is a very straightforward way 
of obtaining the phase diagram, however, to understand 
the physics behind it let us take a step back. When J2 
is relevant, pair-tunnel events in junction 2 will prolifer- 
ate. This means that any field felt by the gas particles of 
type 2 (c;orrc;sponding to the pair-tunnel events in junc- 
tion 2) will be screened by type-2 particles attracted to 
the source of the field. So every type-1 gas particle will 
acquire a screening cloud of particles of typc-2, so that no 
field from the original type-1 particle is felt in junction 
2. To make this a quantitative statement, a type-1 gas 
particle with charge qi exerts the field 



log 



on the type-2 particles. Type-2 particles will then form 
a screening cloud of charge q2 so that: 



which leads to: 



= -Q2El^) = 2g2 log 



92 = Ql . 

R2+r 

now, the field that a test charge of type-1 would feel is: 



n p(ll) 



= -29i^(^?i+^-r-^)log|^| 

log I - 1 , 



(D24) 



H-- 



which we see gives exactly the same result as (D23). In- 
deed this way is more complicated, however it could also 
be employed in more complicated setups, and gives some 
insight as to what physically happens to the system. In 
this case, the charge timncling from lead-1 to the grain, 
partially relaxes through the superconducting junction- 
2. The physical interpretation of the above results is also 
discussed in Sec. (IIIC). 

Next, let us consider the case of relevant J+. Here we 
need to set Ai = — A2 = A. This gives a free energy 
kinetic part: 



S+ 



-I 



Rq 1^ 

27r 277 i?i 



(D25) 



dJi ^ 



1 



1 

Rq 



R1R2 
Ri + R2 



(D26) 



this shifts the phase boundary between the SC* phase 
and FSC phase (Fig. 4(b)) to: 



R1R2 

Ri + R2 



Rq- 



(D27) 



Here too, we can follow the screening principal to get 
the answer. The idea would be that a pair-tunnel event 
would acquire a pair-tunnel couple screening cloud so 
that other pair-tunnel couples won't feel any field. 

This method along with the Self-Consistent-Harmonic- 
Approximation^^ can be used to obtain more insight 
about the behavior of the system. 



APPENDIX E: COULOMB GAS OF 
PHASE-SLIPS REPRESENTATION OF THE 
STRONG COUPLING CASE 

1. Villain transformation - phase-slips 

To treat the strong coupling limit J\, J2 '> 1, we 
need to derive a description of the two-junction system 
in terms of phase-slips: events in which the phase of one 
of a Josephson-junction tunnels from one trough of the 
Josephson cos potential into an adjacent trough. This 
event leads to a voltage drop across the junction (from 
h(j)/2e = V) and hence to dissipation. To derive this 
action, we make use of the villain-transformation. 

Starting with the action (26) 



J dui Rq 



A2 M(l + 4^) 

^1 



27r 27r I 1 fli+fl^+illJla 



lA2 M(l+^) 



2A1A2 



Sj 
(El) 



with: 



Sj 



d- (-JcosAi - JcosAa) (E2) 



Here too we modified the sum over frequencies into an 
integral, and introduced a high-frequency cut-off. 

The assumption of strong J allows us to perform a 
Villain transformation: 



exp (/ drJ (cos(Aj) - 1)) 
Eexp(-/dT| (A, + 27r7?y') 



(E3) 



where r?^^^ maps imaginary time to the integers, and the 
sum on the RHS is over all these functions. The function 
vIt) specifies in which trough of the potential Ji cos Aj 
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the i'th junction is. The essence of the villain transfor- 
mation is that it completely eliminates the dynamics of 
intra-trough motion and only considers the tunneling be- 
tween troughs. The intra-trough near-minimum motion 
is encoded into what will become the fugacity of a phase 
slip, C^- 

It is actually better to use the Fourier transform of the 
time derivative: Fr(r)^^j) = —iuirjlu,) = P\ui)- Incorporat- 
ing this allows us then to write: 



exp {-Sj) « E <^xp - / 



duj J_ 
2tt 2 



+ 



A2 + 27r^ 




Ai -h27r 



—ioj 



(E4) 

Expanding the square and putting it all in the action (26) 
gives: 

c ^ r do. / A 2 / J , l'-l(i+f ) 



R1+R2+ 



2A1A2 



1^1 



R1+R2+ 



+Ai + :li2lPl)l + A2 + :li2lP2}l] . 

(E5) 

Recalling: Z = J2 J -^[Ai]_D[A2] exp(— 5), we are ready 



1, 2 



to integrate out Ai , A2 and get the partition function for 
the phase-slips gas. After doing this and taking the limit 
of large J we get the following partition function: 

Z = E exp / If [i2.p,f ^^^^g^ 



1+- 



l-^l Ri+R2+^^ 

+2(2.pO(2-P2)^^^^^^ 



(E6) 



This is a partition function for a gas that consists of two 
kinds of particles, with pi, p2 being the densities of the 
two gasses. When carrying out the co integrals wc get the 
interaction energy between the gas particles. They are 
(£;('J) is the energy of interaction between two positive 
particles, one of species i and the other from species j): 



-' r..R2\^ H^)- 



E 



(12) 



(E7) 

As in the weak coupling limit, here too we derived a 
Coulomb gas description of the action (26). However in 
the strong coupling limit the gas particles are phase slips, 
which produce a voltage drop over the junction. 



2. From the Coulomb gas to sine-Gordon 

The interacting gas of phase sHps described in (E 1) can 
be encoded into a new sine-Gordon theory, conjugate to 
the original theory (26). It is given by: 

z = jD[e,]jD[e2] 

{~It2^ii^ + Ri)ol + i^ + R2) el - 2re,e2) 

+ / (Ci cos(0i) + C2 cos(^^2) + C- COS (01 - 02))) 

= jD[e,]jD[e2] 

exp (- / ^e^Ge + J ^ (Ci cos(0i) + C2 cos(02))) 

(E8) 

where ^1,2 play the role of fugacities of the phase-slips 
on junctions 1 and 2. By expanding this Sine-Gordon 
theory in the ('s and following the steps of Appendix 
D 1 we recover the Coulomb gas described in (E7). 
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